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Phase change material slurries (PCS) can be employed in the refrigeration and air conditioning systems 
as both secondary refrigerant and cold energy storage media simultaneously, which benefits not only 
the system efficiency improvement by the high cold carry capacity but also the reduction of the 
environment-negative-impact gas emission. This paper reviews the previous researches and develop- 
ments on the applications of PCS, including ice slurry, microencapsulated phase change material slurry 
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(MPCS) as well as clathrate hydrate slurry (CHS), to the secondary loop refrigeration and air 
conditioning systems. The paper mainly focuses on the generation approaches and storage strategies 
of these slurries which are considered as the most important issues relevant to the application of PCS, 
while the flow and heat transfer characteristics of PCS in both tubes and heat exchangers are also 


ae summarized. Moreover, several application cases of PCS and the corresponding operation performances 
CHS are presented. 
© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction 


The application of the secondary refrigerant to the refrigeration 
and air conditioning systems can reduce the charging amount of 
the environment-negative-impact refrigerant, like CFCs and HCFCs, 
etc., consequently the emission of these refrigerants to the atmo- 
sphere is reduced. Meanwhile, the cold energy storage is often 
employed in the refrigeration and air conditioning systems to shift 
the peak-load to the off-peak time, and thus many merits can be 
achieved: reducing the system scale, minimizing the partial load 
compressor operation to improve its efficiency, improving the 
system COP by the extended night time operation at lower 
environment temperature, reducing the financial cost by producing 
and storing the cold energy during night time when the electricity 
is cheaper and releasing the cold energy during day time when the 
electricity is more expensive. The inspiring aspect is that both the 
roles of the secondary refrigerant and cold energy storage medium 
can be undertaken by phase change material slurries (PCS) simul- 
taneously, such as ice slurry, microencapsulated phase change 
material slurry (MPCS) and clathrate hydrate slurry (CHS). The 
latent heat of the fusion of the solid particles containing in PCS 
makes it possible to preserve cold energy with high cold storage 
density, therefore the volumes or sizes of the utilized storage tank, 
piping, pump and the refrigerator are relatively small; the easily 
pumpable performance of PCS through tubes as well as heat 
exchangers without clogging provides the qualification for PCS 
being a good secondary refrigerant with high cold carry capacity, 
then the pumping power is significantly reduced since the required 
flow rate is low. Therefore, PCS was considered promising as 
“functional thermal fluids” [1]. 

The utilization of ice slurry has been widely investigated over 
the world from the fundamental thermal properties to the 
practical applications [2-4]. Ice slurry is a homogeneous mixture 
of very small ice crystals, usually 0.1-2 mm in characteristic size, 
and the base liquid can be either pure water or binary solution of 
water and freezing point depressing agent, e.g. ethanol, glycol or 
salt, etc. The latent heat of ice crystals is about 333.55 kJ kg—', 
which results in a high cold carry capacity. Compared to single- 
phase water or glycol solution, ice slurry with 30 wt% ice fraction 
offers an 80% reduction of flow rate and therefore the pumping 
power while maintaining the same cooling load [5]. Due to its 
excellent cooling and transporting features, ice slurry can be 
employed in many fields, such as the building cooling, super- 
market display, food industry, brewery, mine cooling, fish pro- 
cessing, as well as medical cooling. However, the practical 
application is still limited by the following disadvantages: the 
evaporation temperature is low for air conditioning using; heavy 
mechanical power is required to generate small and smooth ice 
crystals and to maintain the ice slurry in homogeneous state; the 
pressure drop and pumping power of ice slurry with high ice 
fraction are very high. Therefore, the utilization of ice slurry is not 


common except in Japan, and only few projects exist in USA and 
Europe up to date. 

MPCS is a kind of solid—liquid suspension consisting of the carrier 
liquid and small particles of phase change material (PCM) with a 
thin shell. The diameter of the microcapsules is in the range of 
1-1000 um, the core PCM can be solid, liquid, or even mixture. 
Paraffin waxes are often used as the core materials due to its 
extended melting temperature range. The shell material is usually 
natural or synthesized polymer, such as polyester and polyethylene, 
which has high strength and flexibility. The carrier liquid is usually 
water. Compared to non-encapsulated PCM slurry, MPCS owns 
several advantages: the agglomeration of PCM particles and the in 
turn clogging in piping are avoided; the PCM does not contact with 
environment directly, which avoids disadvantageous changes of 
PCM composition and property; the microencapsulation offers a 
protective structure element which improves the mechanical stabi- 
lity, thus permits the frequent changes in volume of the core 
material as phase change occurs; large heat transfer area is achieved 
by the surface of the microencapsulation. So far, however, the 
investigations about MPCS are mostly on its fabrication methods 
and thermal properties [6-8] and the flow and heat transfer 
characteristics through tubes [9-15]. Due to the high cost of MPCS 
preparation, only few laboratorial applications of MPCS to chilled 
ceiling panels [16,17] and an industrial cooling project at Narita 
Airport in Japan [18] are reported in the open literature. 

Clathrate hydrate can be formed in some aqueous solution of 
tetraalkylammonium salts with simple anions (such as halides, 
sulfate, formate, etc.) [19] at the atmospheric condition, in which 
the water molecules form structure lattices as host and molecules 
of another substance fill the lattices as guest; also, gases, like CO2, 
N2, CFCs and so on, can form hydrate in water under the higher 
pressure and low temperature conditions; however, such type of 
gas hydrate slurry would not be discussed in the present paper 
because its application require more rigorous conditions. 
Researchers in Japan [20] applied CHS to the air conditioning 
system as cold storage and transportation medium, and found 
that tetrabutylammonium bromide (TBAB) CHS was quite suita- 
ble to the air conditioning application due to its following merits: 
TBAB CHS can be easily generated at atmospheric pressure and 
room temperature by cooling down the aqueous solution; TBAB 
CHS has a freezing point between 0 and 12 °C depending on TBAB 
concentration in the aqueous solution, which fits quite well with 
the necessary coolant temperature in the air conditioning system; 
the cooling capacity of TBAB CHS is about 2-4 times as large as 
that of chilled water within the temperature range of 5-12 °C; the 
fluidity of TBAB CHS is also good. Researchers in JFE Engineering 
Corporation also successfully applied TBAB CHS to a 4000 m° 
office building [21] for cooling and the operational results showed 
47% reduction of the pumping power. Except the application 
projects conducted by JFE Engineering Corporation [21,22]; however, 
other studies of TBAB CHS still focused on the fundamental 
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Nomenclature 


G specific heat (kJ kg~! K7!) 

D diameter (m) 

dp particle size (m) 

f friction factor (dimensionless) 

fr Fanning friction factor (dimensionless) 

Gz Graetz number (dimensionless), Gz=RePrDy/x 
h heat transfer coefficient (W m~? K7!) 

AH enthalpy change (kJ kg~!) 


H height (m) 

K’ fluid consistency coefficient (dimensionless) 
L length (m) 

n flow behavior index (dimensionless) 

Nu Nusselt number (dimensionless) 

Pr Prandtl number (dimensionless) 

AP pressure drop (Pa) 

dw wall heat flux (W) 

Re Reynolds number (dimensionless) 

Reve modified Reynolds number (dimensionless) 
Ste Stefan number (dimensionless) 

T temperature (°C) 

AT temperature difference (°C) 

v flow velocity (m s7!) 

Xs particle fraction (wt%) 


H viscosity (mPa s) 

A thermal conductivity (W m7! K7!) 
p, ® volume fraction (dimensionless) 
Subscripts 

a average value 

B Bingham 

BT Bingham parameters in Tomita equation 
cf carrier fluid 

f fluid 

h hydraulic 

i inner 

is ice slurry 

local local value 

m mean value 

overall overall value 

S slurry 

w wall 


Ax, particle fraction difference (wt%) 
x distance from inlet (m) 


Greeks symbols 


properties [23], the flow and heat transfer characteristics of TBAB 
CHS flowing through tubes [24-31] as well as in heat exchangers 
[32,33]. In addition, other CHS materials, such as trimethylo- 
lethane (TME) hydrate slurry [34-37] and tetrabutylammonium 
fluoride (TBAF) hydrate slurry [38], have also been studied as the 
thermal energy storage medium. 

This paper presents a review of the current main applications 
of PCS, including the aforementioned ice slurry, MPCS and TBAB 
CHS, to the refrigeration and air conditioning systems as both the 
cold storage medium and secondary refrigerant. The generation 
methods and storage strategies which are very important to PCS 
application are summarized, and the flow and heat transfer 
characteristics of PCS in tubes and heat exchangers are also 
discussed. Some application cases of the practical projects or 
laboratorial studies are displayed, and the corresponding opera- 
tion results are also discussed. 


2. Generation approaches and storage strategies 
2.1. Generation approaches of ice slurry and TBAB CHS 


The generation processes of ice crystals in water or aqueous 
solution and that of TBAB clathrate hydrate are similar in that the 
crystals initiate and grow up in the solution. The generation 
process generally include the following steps: supercooling of the 
base liquid, nucleation and growth of the crystals. The base liquid 
is first cooled down to the temperature lower than its phase 
change temperature but without phase change, the state of which 
is named as the supercooled state. Then the nucleation occurs, 
which is usually heterogeneous since the alien dust particles as 
well as the rough container wall provide nuclei for the nucleation. 
And then it is followed by the growth of crystals. After these three 
steps, size and shape of the crystals are affected by the processes of 
attrition, agglomeration and ripening [39]. If the generated slurry is 
agitated or collides with other materials, attrition will happen due 
to shear stress on the crystals, and large crystals will break off 


during the attrition. Agglomeration may occur during the slurry 
generation or storage without stirring, and small crystals collide 
and stick to each other to form the large size crystals. If the slurry is 
stored for long time, small crystals may dissolve slowly and then 
the size distribution changes, which is called Ostwald ripening. 

Usually, small and globular crystals are expected for easy 
pumping through tubes as well as heat exchangers, therefore 
the attrition process is commonly utilized for crystal size reduc- 
tion and edge smoothing. On the other side, it was reported that 
the ice crystals produced in pure water were dendritic and had 
very rough surface, which resulted in the poor fluidity and easy 
agglomeration and clogging even though the crystal size was 
small [4]. The improvement on this dendritic structure can be 
achieved by smoothing the crystals in chemical way through 
adding the freezing point depressing agents to water, such as 
ethanol, glycol, propylene glycol and sodium chloride and so on. 
Hayashi and Kasza [40] produced ice slurry in glycol solution, and 
found that even very small amount of glycol could improve the 
fluidity of the generated ice slurry significantly, and the higher 
concentration the better fluidity. 

The generation approaches are the pivotal technology for the 
applications of both ice slurry and CHS, since not only the crystal 
properties, i.e., size, shape and roughness, are determined by the 
generation process, but also the system energy efficiency and 
reliability mainly lie in the generator performance. The genera- 
tion of ice slurry has been widely investigated in the alternative 
ways, which can be mainly categorized to the methods using 
moving parts and the methods without moving parts. Therefore, 
the following contexts are mainly focused on the research results 
of ice slurry, and these approaches can also be applied to CHS 
because of the similar generation mechanism as aforementioned. 


2.1.1. Scraped surface generator 

Scraped surface ice slurry generator has been widely investi- 
gated and developed over the world, and it has been applied in 
large-scale systems [5]. The scraped surface generator commonly 


5024 P. Zhang, Z.W. Ma / Renewable and Sustainable Energy Reviews 16 (2012) 5021-5058 


utilizes a shell-and-tube heat exchanger in which the evaporating 
refrigerant flows through the shell side while the crystals are 
generated in the tube side. Driven by a motor, a rotatable scraper 
is mounted on the center of the inner tube to prevent the crystals 
from sticking to the cold heat transfer surface, as shown in 
Fig. 1(a) [39]. Stamatiou [41] and Stamatiou and Kawaji [42] 
developed a compact heat exchanger with scraper inside the 
channel to generate ice slurry, as shown in Fig. 1(b). The scraper 
moves up and down to eliminate the ice crystals on the tube wall. 
The authors obtained 30% power saving by using compact heat 
exchanger instead of tubular one, since four times lower agitation 
speed was required to maintain the ice-free heat transfer surface. 

The blades of the scraper may be plastic or brush type so that 
they can have good contact with the tube wall to detach the 
crystal layer, since it is believed that the nucleation took place on 
the coldest area inside the tube, i.e., the inner tube wall, and then 
the crystals layer grows on the wall forming a thermal resistance 
layer which reduces the heat transfer coefficient between the 
refrigerant and fluid inside the tube. However, other researchers 
thought that the crystals were formed in the bulk of the super- 
cooled liquid but had a tendency to adhere to the cold wall by the 
cohesion force, and this was why the blades were needed. More- 
over, the scraper movement also contributes to the great enlarge- 
ment of the fluid turbulent degree, therefore the heat transfer 
performance of the heat exchanger is promoted. Rodriguez 
Pascual et al. [43] investigated the influences of the plastic blade 
geometries on crystal layer removal efficiency and the slurry 
mixing effect by visualization and simulation. The authors found 
that the mixing of scraped ice crystals with the bulk solution by 


b o 


Staggered blades Chilled surfaces 


move up and down AIN 
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Ice-slury channel 
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Fig. 1. Schematic diagram of the scraped surface ice slurry generators, (a) shell- 
and-tube type [39]; (b) compact heat exchanger type [41] and (c) auger type 
scraper [47]. 


the streamlined blade was inadequate, while the vertical blade 
promoted the mixing but accumulated the crystals in front of it. 

Alternatively, Lakhdar et al. [44] employed stainless steel blades 
rather than the plastic ones. The authors permitted a small gap 
between the blade and the tube wall, and found the permanent ice 
layer only existed when the gap was large enough, about 1 mm 
according to their study. This fact may indicate that the ice crystals 
nucleate inside the solution. Helical-blade auger can also be 
utilized as the scraper for ice slurry generation [45-47], as shown 
in Fig. 1(c). In addition to the effect of scraping the heat transfer 
surface, the rotating auger also provides a negative pumping 
pressure to suck the water (or solution) flowing through the 
generator, hence it can provide the function of the pump. 

The generation of ice slurry using scraped surface is considered 
as very reliable, even in the case of little freezing point depressing 
agent being added to water. During the generation, the increase of 
ice fraction is usually implemented by adjusting the liquid flow 
rate. Using a Sunwell scraped surface ice slurry generator, Wang 
and Kusumoto [5] produced 0-35 wt% ice slurry with about 
100 um ice crystals, and the used base liquid is an aqueous 
solution with 6-10% ethylene glycol. However, drawbacks exist 
in this generator: the tight contact of the plastic scraper blade on 
the tube wall causes strong friction, and then the blades need to 
be replaced after a certain time interval of utilization; since the 
generator employs additional mechanical components, the initial 
investment and the energy consumption are high. 

The scraped surface generator is also suitable to TBAB CHS 
generation. Douzet et al. [48] successfully applied a scraped 
surface heat exchanger (about 200 L, supplied by HeatCraft Ltd., 
France) to TBAB CHS air conditioning system. Daitoku and Utaka 
[49] investigated the scrape effect of different scraper shapes on 
the detachment of TBAB clathrate hydrate from heat transfer 
surface. The authors found that type B TBAB clathrate hydrate 
cohered to the surface more tightly and needed stronger force to 
detach, and a relief angle of the scraper was required to avoid the 
accumulation of crystals on the clearance between the scraper 
and the heat transfer surface when detaching type B hydrate. 


2.1.2. Orbital rod ice slurry generator 

Another type of ice slurry generator using moving parts is the 
orbital rod generator, as shown in Fig. 2 [50]. This generator is 
usually a vertical shell-and-tube heat exchanger with evaporating 
refrigerant flowing from bottom to top in the shell side, while inside 
the shell is a tube bundle consisted many thin tubes. As shown in 
the figure, each of tube has a rotatable whip rod inside, which does 
not rotate in the tube center but around the tube circumference. The 
rotating rods undertake several roles: pushing part of the upper 
income liquid flowing around the tube circumference to form a 
falling film; enhancing the film coefficient and the flow turbulence 
and therefore the heat transfer coefficient between the film and the 
refrigerant; preventing the generated ice crystals from sticking to 
the heat transfer surface, just like the function of scraper in the 
scraped surface ice slurry generator. The rotation of the whip rod is 
much faster (e.g., 850 rpm) than that of the scraper (e.g., 450 rpm), 
however, it is believed that the falling film acts as the lubricant 
which minimizes the drive power of the motor. 

Orbital rod ice slurry generator is also widely used, e.g., Paul 
Mueller Company [51] developed ice slurry units using this 
generator which has the standard production sizes from 3 to 
520t (about 10.55-1828.68 kW). They utilized the generator in 
many fields, including refrigeration, HVAC and process cooling and 
so on. Stamatiou et al. [39] compared the performances of the 
orbital rod and the scraped surface ice slurry generators. It is 
believed that both the two kinds of generators are stable, and the 
generated ice crystals are homogeneous due to influence of the 
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Fig. 2. Schematic diagram of the orbital rod ice slurry generator [50]. 


rapid rotated scraper or rods. The scraped surface ice slurry 
generator has a slight higher cooling capacity, however, the drive 
power for the scraper rotation is also larger; the overall heat 
transfer coefficient of the orbital rod ice slurry generator is higher 
due to the falling film heat exchange. The whip rod is metallic and 
has no direct contact with the tube, therefore it survives longer 
time than that of the plastic scraper blade which attaches to the 
tube tightly as aforementioned. However, the drive mechanism of 
whip rods is complicated and there are many associated problems. 


2.1.3. Fluidized bed generator 

The fluidized bed heat exchanger was proposed in 1970s for 
sea water desalination [52]. Small solid particles are inserted into 
the tubes in a shell-and-tube heat exchanger to detach the 
adherence on the tube wall. The insert particles are usually 
cylindrically shaped stainless steel made by wire, or globular 
glass balls. As the fluid flowing upward in the tube, the solid 
particles are fluidized and continuously impact on the heat 
transfer wall, therefore the deposits on the wall are broken off 
and fall into the fluid, as shown in Fig. 3 [53]. Meanwhile, the heat 
transfer coefficient is also improved by the following three 
mechanisms [53]: the laminar boundary layer that acts as heat 
transfer resistance layer is broken up by the fluidized particles; 
the solid particles transport the heat from the heat transfer wall 
to the fluid by conduction as the particles contacting the wall and 
then the fluid; the deposit layer on the wall is broken off. 
Therefore, the heat transfer coefficient is determined by the flow 
velocity of the fluid and also the bed voidage (the ratio of the 
space that not occupied by the solid particles to the entire space). 

Klaren and van der Meer [54] first applied the fluidized bed 
heat exchanger to produce ice slurry in 1991, then Meewisse and 
Infante Ferreira [55,56] and Pronk et al. [53,57-59] conducted a 
lot of works on producing ice slurry by the fluidized bed. 
Schematic diagrams of the typical fluidized bed ice slurry gen- 
erators are shown in Fig. 4. The generator shown in Fig. 4(a) is the 
stationary fluidized bed, in which the insert solid particles are 
kept inside the tube during the operation. Shown in Fig. 4(b) is the 
so-called circulating fluidized bed, where the solid particles will 
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Fig. 3. Schematic diagram of the effect of fluidized solid particles [53]. 


flow out of the tube together with the generated ice slurry, and 
then are separated from the slurry and return to the bottom via a 
downcomer. Measurement results by Meewisse and Infante 
Ferreira [56] showed that the heat transfer coefficients of the 
fluidized bed during ice slurry generation was in the range of 
2500-4000 W m~? K~!, which was about 2-4 times as large as 
that of single-phase flow and was not influenced by the ice slurry 
generation. They also proposed the empirical heat transfer corre- 
lations which predicted the heat transfer accurately. Compared to 
the stationary mode, the use of the circulating fluidized bed is 
more flexible since the solid particle can flow at a high velocity 
range, and the advantages include the improvement of the ice 
crystals breaking off from the tube wall by increasing flow 
velocity and also the heat transfer coefficient is increased by 
1-25% [58]. 
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Fig. 4. Schematic diagram of the fluidized bed ice slurry generator, (a) stationary type and (b) circulating type [57]. 


The shell-and-tube design of the fluidized bed is simple and 
results in an inexpensive and low operation cost for ice genera- 
tion. Due to the violent impacts of solid particles on the tube wall, 
the generated ice particles can be homogeneous, which is usually 
in the shape of cylindrical or ellipsoidal with the size of about 
300 um [57]. The main disadvantage of the fluidized bed ice 
slurry generator is the existence of a maximum allowable tem- 
perature difference between the fluid inside the tube and the 
refrigerant (or the tube wall), below which ice slurry can be 
generated stably with nearly constant heat transfer coefficient. 
The higher temperature difference enlarges the amount of ice 
crystals sticking on the tube wall, which are hard to be removed 
efficiently by the impact of solid particles. The gradually thicken- 
ing crystal layer deteriorates the heat transfer between the slurry 
and the refrigerant, thereafter ice slurry generation becomes 
unstable. Meewisse and Infante Ferreira [55] measured this 
temperature difference limit in a circulating fluidized bed, as 
shown in Fig. 5. It was found that the temperature difference limit 
varied linearly with the freezing temperature of the aqueous 
solution. Lower freezing temperature (or higher additive concen- 
tration) caused a higher temperature difference limit. The other 
disadvantage of this ice slurry generator is the requirement of 
higher fluid velocity, so that the momentum s of the solid particles 
are sufficient to impact the wall for ice crystal layer detachment. 


2.1.4. Supercooling method 

Ice slurry generation using supercooling concept employs the 
following steps: water or aqueous solution is cooled in an 
evaporator to a few degrees below its freezing temperature 
without crystallization; the supercooled liquid is forced to flow 
through a special device, i.e., supercooling releaser, to release the 
supercooling state; then ice crystals are generated in the liquid 
and the liquid turns to solid-liquid ice slurry. This method is 
considered as simple and having low operation cost since there is 
no other complex apparatus. However, the evaporator for super- 
cooling and the supercooling releaser must be carefully designed 
to achieve stable ice slurry generation. 

The heat exchanger for supercooling must ensure that there is 
no crystallization happened inside, otherwise the slurry genera- 
tion may stop due to the blockage by the accumulatively gener- 
ated ice crystals. Bédécarrats et al. [60] regarded that the 
stochastic character of crystallization resulted in a hard control 
of ice generation within heat exchanger and it was necessary to 
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Fig. 5. Range of stable ice slurry production using circulating fluidized bed, 
(a) range of NaCl and (b) comparison of the operating ranges of NaCl and ethylene 
glycol (EG). (ATmax is the temperature difference between the slurry and the 
refrigerant, the bed voidage is 0.8) [55]. 


control the evaporation temperature accurately. They studied the 
influences of the supercooling degree, refrigerant temperature 
and fluid flow rate on the ice generation. By developing models 
for predicting the fluid temperature variation, they tried to avoid 
the ice generation inside supercooling heat exchanger. 
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Several requirements for the supercooling heat exchanger are 
given as follows: 


(1) The heat transfer surface should be smooth, since the rough 
surface enhances the crystallization. Through the observations of 
crystallization within tubes with different roughness (between 
0.63 and 13.3 um), Faucheux et al. [61] found strong influence of 
surface roughness on supercooling degree and concluded that 
larger roughness caused lower supercooling degree. 

(2) The flow channel for the water or solution should be regular, 
since the sudden change of the channel, e.g. bend, contraction 
or expansion, results in strong collision between liquid and 
solid surface, which accelerates the nucleation. 

(3) The heat exchange inside the evaporator should be uniform 
everywhere of the heat exchanger; otherwise crystallization 
may occur locally where the temperature of the fluid is low. 


The supercooling releaser is also important to the effective 
crystal generation, and several releasing approaches are investi- 
gated and utilized as follows: 


(1) The collision of the supercooled fluid with a solid surface, or 
free liquid surface or liquid bulk in the tank, or with itself, as 
shown in Fig. 6 [62-65]. The collision provides the momen- 
tum increase to release the supercooling state. The probability 
of freezing mainly depends on the collision force [66]. More- 
over, the system requires special device to avoid ice crystals 
agglomeration near the impact place. 

Ultrasonic vibration: Researchers [67,68] reported that the ultra- 
sonic vibration can strongly promote the transition of super- 
cooled water to ice slurry, and this method is reliable for ice 
generation but it is only limited to laboratorial investigation. 
Partial low temperature: Partial liquid being cooled to the 
temperature much lower than its freezing temperature can be 
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helpful to the nucleation. This process can be achieved by a 
small electronic refrigerator. 

(4) Electric charge: Though scientists have tried to investigate the 
mechanism of the nucleation by electric charge, the basic 
mechanism is still controversial [69]. However, some experi- 
mental results have been identified to shown the influential 
factors. Okawa and Saito [66] disclosed supercooling released by 
the utilization of electric charge with the voltage less than DC 
100 V, and found that the influential factors were the voltage, the 
size of electrode and the distance between the two electrodes. 


Okawa and Saito [66] also treated the reciprocally rubbing of 
glass against glass as one of the effective supercooling release 
approach, which was depended upon the frequency and ampli- 
tude of the rubbing. On the other hand, the authors believed the 
convection, vibration, and non-contacting shock had no effect on 
supercooling release. 

Tanino and Kozawa [70] and Kozawa et al. [71] applied a pipe 
type releaser for large-scale ice slurry generation, which generated 
2.5 wt% ice slurry for one circulation, and they finally obtained 
35 wt% ice slurry. While Ise et al. [72] used this kind of releaser to 
generate about 23 wt% ice slurry for a building usage. 

Nearly the same approach can be employed for CHS generation, 
e.g., JFE Engineering Corporation [73] successfully applied TBAB CHS 
air conditioning to several large-scale buildings [21,22,74]. Ogishi 
et al. [75] employed a tube type heat exchanger (such as shell-and- 
tube heat exchanger) as the supercooling evaporator for TBAB CHS 
generation. The heat exchanger also took the role of supercooling 
releaser, therefore there was at least one of the tubes had irregularity 
for supercooling release. Hydrates would be generated inside the heat 
exchanger. The authors reported the deposition of hydrate crystals on 
the tube inner wall and treated them to be useful nuclei for stable 
hydrate generation, while the other hydrates were removed by the 
shear force of the flowing aqueous solution or slurry. The thickness of 
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Fig. 6. Collision modes for supercooling release, (a) supercooled liquid and solid surface [63]; (b) supercooled liquid and liquid surface in the tank [64]; (c) supercooled 


liquid and bulk liquid in the tank [65] and (d) supercooled liquid and itself [62]. 
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the deposition layer was further controlled by adjusting the slurry 
flow rate and regulating the refrigerant temperature. Moreover, Song 
et al. [76] validated the effect of ultrasonic vibration on TBAB CHS 
generation in 25 wt% TBAB solution, and the nucleation was observed 
to be generated in a very short time under the ultrasonic vibration 
and the crystallization process also finished in several seconds. The 
authors found that the crystals may be cracked under the force of 
ultrasonic vibration, and the obtained TBAB hydrate crystals were 
fine, smooth and uniform without agglomeration. 


2.1.5. Direct contact heat exchanger 

Using the direct contact heat transfer between water or aqueous 
solution with immiscible cool refrigerant, including primary and 
other secondary refrigerants, ice slurry can be generated directly 
inside the storage tank [77-83]. Fig. 7(a) shows the utilization of 
primary refrigerant. During the evaporation of the two-phase primary 
refrigerant, dispersed ice particles are generated while the refrigerant 
can be easily separated since it turns to vapor. However, the heat 
exchanger needs to work under high pressure. Fig. 7(b) shows the 
application of secondary refrigerant which has lower freezing tem- 
perature and usually larger density than water, the refrigerant sinks 
to bottom of the tank after heat exchange and then returns to the 
pump and chiller for another cycle. 
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Fig. 7. Schematic diagrams of the direct contact heat exchanger for ice slurry 
generation, (a) using primary refrigerant [78] and (b) using secondary refrigerant [83]. 


Within a storage tank of 0.3 m in diameter and 1.1 m in height, 
Kiatsiriroat et al. [78] utilized the primary refrigerant, including 
R22 and R12 in the temperature range of —8 to —15 °C, as the 
coolant for direct contact heat exchange with water. The flow rates 
of the two refrigerants were 0.02-0.08 kgs7! (about 0.28- 
1.13 kg m~? s7!) while water was stagnant and kept inside the 
tank. They found that the volumetric heat transfer coefficients in 
the cases of R12 and R22 were around 2-16 kW m~? K`! and 
2-52 kW m~? K7}, respectively. The promoted heat transfer per- 
formance by R22 was attributed to its higher water solubility. 
Thitipatanapong et al. [79] investigated ice generation using R12 
and R134a in direct contact heat exchange while both the refrig- 
erant and water were forced flowing into the heat exchanger. It 
was noticed that higher water flow rate resulted in higher volu- 
metric heat transfer coefficient. The volumetric heat transfer 
coefficients obtained by the authors was about 91.62 kW m~? K7! 
and 209.7 kW m~? K~! in the cases of R12 and R134a with water 
flow rates of about 16.35kgm~?s~' and 16.20 kg m~? s™!, 
respectively, which was higher than those obtained by Kiatsiriroat 
et al. [78]. Due to the saturation pressure of refrigerant, the ice 
slurry forming temperature increased to 4.0°C and 8.8 °C in the 
case of R12 and R134a, respectively, therefore the storage tem- 
perature at about 4.5 °C and 10.0 °C can be achieved, which was 
considered suitable to air conditioning system application. Thong- 
wik et al. [80] injected low temperature CO2 (— 15 to —60 °C) to 
water at a flow rate of 0.003-0.017 kg s~', and the corresponding 
volumetric heat transfer coefficient was ranged from 1 to 
25 kW m~? K~—!. After the direct contact with the primary refrig- 
erant, the base liquid is inevitably contaminated by the compressor 
oil, therefore the authors added the compressor oil and surfactant 
Tween 60 (polyoxyethylene sorbitan monostearate, which was 
used to merge water and oil) to water to investigate its influence 
on slurry generation. It was found lower oil concentration 
increased the risk of blockage around injector since the non- 
merging water formed plain ice and blocked the injector, and the 
optimum composition of water/oil/surfactant was 100/6/1. 

Researchers in Singapore [81-83] proposed a special coolant, 
named Fluroinert (FC-84), as the suitable secondary refrigerant used 
in direct contact ice slurry generation. FC-84 is heavier than water, 
the density of which is about 1730 kgm~3, while the freezing 
temperature is — 95 °C. As shown in Fig. 7(b), FC-84 was sprayed by 
a nozzle at the top of the storage tank. The coolant was at about 
—10 to —15°C while the flow rate was 8-12 L min™!. Through 
modeling and calculation, the overall heat transfer coefficient of 
the heat exchange between water and FC-84 was in the range of 
3.0-6.5 kW m~? K~!. The generated ice crystals were in the shape 
of thin and long flake with the length of about 1-4 mm, which was 
quite large compared to those obtained by other methods. 

Compared to other ice slurry generators, the utilization of 
direct heat exchange avoids the sticking of ice crystals on the heat 
transfer surface and strongly enhances the heat transfer perfor- 
mance. Therefore the COP of the refrigerator increases, for 
example, about 3.4-3.6, higher than that of the conventional 
system [77]. The disadvantages of this method are as follows: the 
nozzle used in the system must be well designed, otherwise 
blockage may occur; small amount of primary refrigerant trapped 
in the ice crystals will be released to the environment. 


2.1.6. Vacuum ice slurry generator 

Ice slurry can be generated in a deep vacuum freezer. Water 
inside the freezer is kept at deep vacuum close to the triple-point 
of water, hence part of the water becomes vapor, which absorbs 
the latent heat of vaporization from the left water, and it is then 
frozen and converted to ice particles or ice-water mixture. 
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Fig. 8. Vacuum ice slurry generator used in thermal energy storage system [85]. 


IDE Technologies Ltd. [84-86] developed a vacuum ice slurry 
generator, and the schematic diagram is shown in Fig. 8. A high 
speed centrifugal compressor driven by a motor was utilized to 
generate deep vacuum condition. The generated ice slurry inside 
the freezer, containing 0.5-1.0 mm ice crystals, was then pumped 
to the storage tank, and the slurry fraction could reach 50%. To 
generate ice slurry continuously, water vapor was compressed 
and then forced flowing through a condenser, where the vapor 
was cooled by chilled water (e.g., 5 °C) and converted to liquid 
phase. The cooling capacity of this generator can achieve 
3500 kW. The generator has been applied to district cooling and 
heating, thermal energy storage, mine cooling and snow making, 
and hundreds units have been installed worldwide [86]. However, 
due to large specific volume of water vapor under deep vacuum 
condition, the size of the generator is quite large. On the other 
side, Bapat and Kulkarni [87] employed a lithium bromide vapor 
absorption machine to generate the vacuum condition for ice 
slurry generation, which is compact in size and less electric 
energy consumption. 

Researchers in Korea [88,89] sprayed water droplets with 
initial temperature at 20°C and 50 um in size by nozzles in a 
vacuum chamber to generate ice slurry. The pressure in the 
chamber was kept lower than the pressure of triple point, and 
spherical ice crystals with size below 300 um were then formed. 
To avoid the generation of ice crystals around nozzles, the authors 
added 7% ethylene glycol to water. Different from the vacuum 
generator shown in Fig. 8, the formed vapor was then removed 
from the chamber. The feeding of water was maintained at low 
rate to obtain water droplets, therefore the cooling capacity was 
only 0.23 kW (theoretical value), which could be improved by 
increasing the nozzle number. 


2.2. Storage strategy of ice slurry and TBAB CHS 


The storage of ice slurry and TBAB CHS in tanks is mainly 
conducted in two modes, homogeneous storage and heteroge- 
neous storage [90], as shown in Fig. 9. 

As shown in Fig. 9(a), the homogeneous storage required 
agitation to maintain the slurry fluidizing inside the storage tank 
and to avoid the agglomeration of particles. The agitation is 
usually provided by propeller, pumps, paddles and so on. Since 
the absence of agglomeration in homogeneous slurry storage, it is 
very easy to extract the slurry with very fine crystals and the 
pumping is also quite reliable. It was reported [91] that the 
maximum ice fraction in homogeneous storage could be up to 
30-35 wt%, beyond which the power consumed by the agitator 
increased significantly, e.g., the power consumed by Danish 
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Fig. 9. Storage strategies of ice slurry, (a) homogeneous storage and 
(b) heterogeneous storage [90]. 


Technological Institute to agitate 52 wt% ice slurry was 10 times 
larger than that of 30 wt%. 

Since the agitation is significant energy consuming, Meili et al. 
[92] reported that the mixing of the slurry could be stopped for 
about 12h, and then homogeneous ice slurry could be still 
obtained after restarting the agitation. Moreover, adding freezing 
point depressing agents to water can smooth the generated ice 
particles and prevent the agglomeration. However, Grandum and 
Nakagomi [93] and Inaba et al. [94] believed the additions of 
freezing point depressing agents resulted in a low refrigerator 
COP because of the depressed freezing temperature. Therefore, 
the authors utilized other additives, for example, antifreeze 
proteins (AFPs), antifreeze glycoproteins (AFPs), silane coupling 
agents, surfactants and so on and considered them as suitable 
since the low additive concentration (less than 1 wt%) led to 
relative high freezing temperature. 

Fig. 9(b) shows the method of heterogeneous storage. Due to 
the different density of the solid crystals from that of the carrier 
liquid, such as ice is lighter than water and TBAB hydrate is 
heavier than the solution, crystals in heterogeneous storage tank 
tend to accumulate at the top or deposit at the bottom of the tank. 
Compared to homogeneous storage, the ice fraction in hetero- 
geneous storage tank can reach about 50-60 wt% [91], and the 
necessary volume of the storage tank is then reduced. The difficult 
thing in this storage mode is that the extraction of ice from the 
tank. Eglof et al. [90] considered using a screw pump to extract 
the high fraction ice slurry from the top of the tank and then 
mixed it with the low fraction ice slurry pumped from the bottom 
region through a special valve. This proposal has not been 
experimentally validated yet and needs further investigation. 
Douzet et al. [48] employed this approach for TBAB CHS storage, 
they pumped the high fraction TBAB CHS (about 40-50% in 
volume) from the bottom of the tank using a lobes pump and 
mixed it with the liquid sucked through the top of the tank using 
another pump, then they obtained TBAB CHS with about 20% 
volume fraction for air conditioning application. Actually, pure 
ice particles can be obtained by a harvester [91], though the 
shaved ice particles (e.g., 1-2 mm) were larger than the original 
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generated ones (e.g., 0.1 mm), the particles are mixed with the 
cold water extracted from the bottom of the tank to form ice 
slurries with different fractions. 


2.3. Fabrication of MPCS 


MPCS requires the technology that coating the PCM with 
suitable shell material to envelop the core PCM. Gibbs et al. [95] 
reported several forms of the microcapsules, as shown in Fig. 10. 
Generally three steps are required for microcapsules preparation: 
forming the shell around the PCM, ensuring there is no leakage, 
and ensuring the undesired materials are not included in the 
capsule. Though the microencapsulation technology was discov- 
ered and developed in the 1940s and 1950s [96] and then many 
approaches were developed for food industrial [97], only coacer- 
vation, polymerization, spray drying were found being employed 
in PCM encapsulation. These processes are simply introduced in 
the following contexts, the detailed information about the fabri- 
cation and the further properties of the formed microencapsu- 
lated PCMs can refer to the review articles [7,8]. 


2.3.1. Coacervation 

There are two types of coacervations: simple coacervation and 
complex coacervation. The complex one is commonly used for 
PCM microencapsulation. First, polymer in water as the coating 
material, e.g. gelatine in water, is prepared. Then the core 
material is dispersed to the gelatine solution at a temperature 
range of 40-60 °C, at which the coating material is liquid. The 
simple coacervation employs only one kind of colloid, while 
another colloid, such as gum Arabic, is added to the solution for 
complex coacervation at this step. The step is followed by 
coacervation. Simple coacervation is induced by change of condi- 
tion, such as the temperature, the addition of nonsolvent or 
microions, which causes the dehydration of macromolecules. 
Complex coacervation is usually induced by increasing the poly- 
mer concentration or reducing the solution pH (usually 4-4.5). 
The system is then cooled to room temperature. Finally, the 
microcapsules are hardened and isolated. Hardening is conducted 
by cross-linking agent, such as formaldehyde, the increasing of pH 
to 9-11 (by NaOH solution, e.g.) and temperature to 50 °C can 
hasten this process. The temperature should be reduced and 
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Fig. 10. Several forms of the microcapsules [95]. 


maintained at about 5-10°C for 2-4h, then the capsules are 
filtered and dried in a vacuum oven. 

Using the complex coacervation, Hawlader et al. [98] obtained 
around 1 um spherical and smooth microcapsules of paraffin 
wax. The latent heat was about 234.05] g~', 224.31Jg~' and 
196.38 J g~' when the core-to-coating ratio was 2:1, 1:1 and 1:2, 
respectively. The encapsulation efficiency was in the range of 
80-95%, and it decreased as increasing the ratio of core-to-coating 
since incomplete coating on core material might occur when the 
coating material was deficient. Ozonur et al. [99] applied the 
same method to encapsulate the natural coco fatty mixture for 
thermal energy storage (melting temperature of the PCM was in 
the range of 22-24 °C). The size of the obtained microcapsules 
varied from 1 um to 1 mm depending on the stirring speed and 
time. After encapsulation, there was an increment of about 7 °C 
on the melting temperature. 


2.3.2. Polymerization 

Polymerization are used widely for fabricating microencapsu- 
lated PCM, and there are several different type polymerization 
approaches. Yang et al. [100] and Brown et al. [101] employed 
the so-called “in-situ” polymerization to encapsulate tetradecane 
and dicyclopentadiene, respectively. The former study utilized 
several types of polymers, e.g., polyvinyl acetate (PVAc), polystyrene 
(PS), polymethylmethacrylate (PMMA), polyethyl methacrylate 
(PEMA). The authors mixed the PCM and polymer with surfac- 
tants, initiator and deionized water in a flask accompanying with 
agitation, and then heated the mixture to and maintained it at a 
certain temperature. The mixture was then cooled down naturally 
by continuous stirring for 4-5 h. It was found that the coating 
material strongly influenced the thermal properties of the formed 
MPCM. After encapsulation, the melting temperature decreased 
from 6.7 °C (pure tetradecane) to 5.97 °C as the coating material 
was PMMA and to 5.68 °C as the coating material was PEMA, 
while the freezing temperature increased from 1.98 °C to 2.95 °C 
and to 3.19 °C, respectively. The obtained microcapsules were 
5-30 um in diameter and 66.26 kJ kg~' in melting heat with 40% 
tetradecane. The latter study described their detailed preparation 
process in Fig. 11. Using this process, 10-1000 um microcapsules 
were obtained by adjusting the agitation rate between 200 and 
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Fig. 11. In-situ polymerization for microencapsulation [101]. 
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2000 rpm, and the wall thickness was in the range of 160- 
220 nm. It was noticed that the surface roughness reduced as 
the agitation rate increased. Fill content of the obtained micro- 
capsules was beneficially up to 83-92 wt%. Using in-situ poly- 
merization and the assistance of ultrasonic, Fang et al. [6] coated 
n-octadecane with polystyrene as shell and obtained very fine 
nanoencapsulated PCM, the particles were regular spheres and 
diameter was only 50-200 nm. The phase change temperature 
was found very close to that of pure n-octadecane, while the 
latent heat dropped nearly 50% due to the existence of the coating 
material. 

Alkan et al. [102] prepared the microcapsules using emulsion 
polymerization. They used docosane as core and PMMA as shell. 
Docosane was first mixed with deionized water and Triton X-100 
(surfactant), and then was heated to the temperature higher than 
the melting temperature of docosane under nitrogen atmosphere. 
This process was followed by the addition of methylmethacrylate 
and allyl methyl acrylate and FeSO,-7H20 solution, of which the 
allyl methyl acrylate was used as cross-linking agent to prevent 
the formation of secondary particles. The mixture was then 
agitated violently for 30 min. Na2S207 solution and tertbutylhy- 
droperoxide solution were added and the mixture was heated to 
90 °C. Thereafter the polymerization started. After casting water, 
the precipitate was dried under vacuum at 40°C for 24h. The 
average diameter of the obtained PMMA microcapsules was 
0.16 um. The influences of microencapsulation on the melting 
and freezing temperature were slight, and only 0.6 °C decrease of 
the melting temperature and 0.2 °C increase of the freezing 
temperature were observed after the microencapsulation. The 
melting and freezing heat were measured, which were 
54.6 kJ kg~! and 48.7 kJ kg~!, respectively. Therefore the slurry 
contained 28 wt% docosane. 

Cho et al. [103] and Siddhan et al. [104] applied almost the 
same interfacial polymerization to encapsulate the PCM (octade- 
cane) using TDI (toluene-2, 4-diisocyanate) and DETA (diethyle- 
netriamine) as coating material. Surfactant NP-10 [poly(ethylene 
glycol) nonylphenyl ether] was treated as emulsifier by the 
former study while the latter study used NP-9.5. The NP-10 (or 
NP-9.5) aqueous solution and organic solution of octadecane, 
cyclohexane and TDI were prepared and mixed with each other 
accompanying with violent agitation to form an oil-in-water 
emulsion. Then the DETA aqueous solution was added to the 
emulsion drop by drop to start the interfacial polymerization 
which happened between TDI and DETA at the oil-water inter- 
face. The rotation speed of the agitator was slowed down and the 
mixture was heated to 60°C. This state was maintained for 
45 min (the latter study added more DETA solution for further 
diffusion and polymerization). Finally, the microcapsules were 
filtered and washed and dried. Both studies obtained small 
capsules, 0.1-1 um in the former study and less than 14 um in 
the latter one. By additionally adding DETA solution, Siddhan 
et al. [104] obtained higher core fill content (about 40-68 wt%) 
than that obtained by Cho et al. (26.0-46.4 wt%). The encapsula- 
tion efficiency was 75.4-87.2% in the former study, while the 
latter study noticed that the efficiency was largely depended on 
the core-to-monomer ratio (varied from about 50% to 90% as the 
ratio increased from about 2 to 6.5) and PCM-to-cyclohexane 
ratio (varied from about 87% to 97% as the ratio increased from 
about 1 to 12). 

Ai et al. [105] used suspension polymerization for n-hexade- 
cane microencapsulation, while polystyrene was used as shell. 
Considering casein in solution and dispersed form as an excellent 
emulsifier, the authors utilized casein instead of the necessary 
surfactants as emulsifier in other polymerization processes to 
emulsify the mixture of core material, monomers and initiators. 
The usage of casein provided extra effect of protecting the oil 


droplets against recoalescence before the microcapsules forma- 
tion. Thus, the application of casein in polymerization process 
could simplify the process since no more additives were needed. 
The size of the obtained microcapsules was 3-15 um. The latent 
heat was in the range of 76.01-124.7 kJ kg~! depending on the 
polymerization time, which was much lower than that of the pure 
material, 236 kJ kg~!. Sanchez et al. [106] also considered this 
method as simple, cheap and technically feasible. They coated 
several kinds of phase change materials, including paraffin wax, 
tetradecane, nonadecane, etc., by a coating material of styrene. 
The other necessary materials included benzoyl peroxide as 
initiator, polyvinylpyrrolodone as stabilizer and water. After the 
violent agitation, the mixture was heated to 110 °C at one atmo- 
sphere for 6h, then it was washed, filtrated and dried. The 
capsules size ranged between 20um and 700um with the 
average value of 237 um, which was also dependent on different 
core materials. The core fill content for paraffin wax, tetradecane 
and nonadecane was 20.56 wt%, 39.10 wt% and 49.36 wt%, 
respectively. 


2.3.3. Spray drying 

Spray drying is a technology that can generate powder 
instantaneously by atomizing the liquid production in a hot air 
or nitrogen [107], which is considered as economic for micro- 
encapsulation. The spray drying process in microencapsulation 
involves the following basic steps: mixing the core solution and 
the coating material solution; heating and agitating the mixture 
to obtain homogeneous mixture; if needed, adding the emulsifier 
to the mixture, which depends on the properties of the coating 
material; ensuring the formed emulsion are stable over a certain 
period; atomizing the emulsion in the drying chamber where the 
solvent (usually water) is evaporated; cooling and dehydrating 
the obtained particles. Hawlader et al. [98] employed a spray 
dryer to spray dry the emulsions of paraffin wax and the coating 
material with a feeding rate of 20 mLmin~! by a centrifugal 
atomizer operated at 25,000rpm. The authors compared the 
microcapsules produced by spray drying with that by complex 
coacervation, and it was found the spray drying made much 
smaller particles (around 0.1 um compared to around 1 um), 
however, the latent heat was also smaller than that produced 
by the complex coacervation (about 221.5 kJ kg~' compared to 
234 kJ kg~!, when the core-to-coating ratio was 2:1). 


3. Flow and heat transfer characteristics 


Flow and heat transfer characteristics of PCS in tubes and heat 
exchangers are significant to its application. The rheology, flow 
and heat transfer behaviors of ice slurry have been previously 
reviewed by Ayel et al. [108] in 2003, and Zhang et al. [7] also 
presented the flow and heat transfer characteristics of MPCS and 
CHS in 2010. However, there are many new researches about the 
PCS that have been conducted thereafter. Therefore we emphasize 
the present review on the recent research details of the flow 
pattern, flow frictions and heat transfer coefficients of PCS in 
tubes and heat exchangers and also the impact of particle fraction 
and tube diameter on mentioned behaviors. 


3.1. Flow pattern of ice slurries 


The usual flow patterns of solid-liquid two-phase flow are 
homogeneous flow, heterogeneous flow, moving bed and station- 
ary bed. The pressure drop increases with flow velocity at the 
homogeneous and heterogeneous flows but decreases at moving 
bed and stationary bed flows. All these flow patterns have ever 
been observed in ice slurry flow, however, Kitanovski et al. [109] 
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pointed that the flow pattern classification by simple observation 
showed quite large uncertainty. 

Using the Ultrasonic Doppler Profiler (UVP), Vuarnoz et al. [110] 
obtained both the vertical and horizontal velocity profiles of 
7.5-7.6% ice slurry flowing through a 20 mm x 20 mm rectangular 
channel, as shown in Fig. 12(a). With the decrease of flow rate, 
both velocity profiles became bulgy in the axes. The vertical 
profiles were more sensitive to the flow rate reduction due to the 
gravity effect, therefore the flow became slower at the top and 
faster at the bottom since the accumulation of ice particles at the 
top, as shown in the figure; meanwhile, the horizontal profiles 
remain symmetrical as decreasing flow rate. Using a hot film 
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anemometer (HFA) and an on-line ice slurry sampling/calorimetry 
technique, Stamatiou and Kawaji [111,112] obtained the velocity 
profiles of 0-5.9 vol% ice slurry and 3.1-9.0 vol% ice slurry with 
heat exchange flowing through a vertical 610 mm x 310 mm x 
25.4 mm channel as shown in Fig. 12(b1)-(b4). When ice fraction 
was low, i.e. 0-0.9 vol% in Fig. 12(b1), the velocity profiles were 
hardly influenced by the ice particles, therefore the results were 
similar to that of single-phase fluid; ice particles influence the 
velocity profiles when the ice fraction is in the range of 0-5.9 vol%, 
as shown in Fig. 12(b2), larger ice fraction results in blunter 
profiles; when heat flux was imposed to ice slurry flow, the 
formation of ice-free layer near the heated wall led to larger fluid 
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Fig. 12. Velocity profiles of ice slurry flow, (a) xs=7.5-7.6%, rectangular channel 20 mm x 20 mm [110]; (b1)-(b4) base liquid was 6.2 wt% NaCl brine, dm=0.1-0.2 mm, 
vertical rectangular channel 610 mm x 310 mm x 25.4 mm [111,112], (b1) adiabatic flow at very low ice fraction; (b2) adiabatic flow for various ice fractions; (b3) diabatic 
flow for various heating powers; (b4) diabatic flow for various ice fractions. The reference and equations in (b1)-(b4) are those in [111,112]. 
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density in this area than fluid in the flow center, therefore the flow 
near the heated wall decelerated as shown in Fig. 12(b3) while the 
flow in the tube center changed little; Fig. 12(b4) shows that there 
is no clear influence of ice fraction on the velocity profiles of ice 
slurry under heating. 

Ice fraction profile of ice slurry is related to many factors, e.g., 
density of carrier liquid, density of ice particle, shape and size of 
ice particle, tube diameter, and the flow velocity [109]. Kitanovski 
and Poredos [113] modeled and numerically calculated the 
concentration profiles of 10 wt% and 25 wt% ice slurry flow 
though 27.2 mm tube, as shown in Fig. 13. It was believed that 
the 10 wt% ice slurry flow was homogeneous at v,=2.0 m s~! and 
was heterogeneous at v,;=0.5-1.5ms~' and moving bed 
appeared at v,=0.5 ms~!, as shown in Fig. 13(a), and the larger 
ice faction resulted in earlier transition from homogeneous flow 
to heterogeneous flow as comparing Fig. 13(a) and (b). 

With the reduction of flow velocity, flow of ordinary solid- 
liquid two-phase slurry will transit from heterogeneous flow to 
moving bed at a critical velocity, which is named as deposition 
velocity. On the basis of velocity profiles, Kitanovski [114] 
evaluated the deposition velocities of ice slurry with different 
ice particle diameters as well as different tube diameters, as 
shown in Fig. 14(a). As the ice fraction increasing from 5 wt% to 
25 wt%, the deposition velocity of ice slurry had a maximum value 
at x,=15 wt%. Both the increase of tube diameter and ice particle 
diameter caused easier deposition corresponding to the larger 
deposition velocity as shown in the figure. Rached et al. [115] 
reported a new method to determine the deposition velocity of 
ice slurry. They measured the cross-section areas of the moving 
bed as increasing the flow velocity gradually, and then extra- 
polated the deposition velocity corresponding to the disappear- 
ance of the moving bed which could not been observed directly. 
Deposition velocities of 10-32 wt% ice slurry flowing through 
50 mm and 75 mm tubes were obtained, as shown in Fig. 14(b). 
As seen from the figure, the deposition velocity reduced with ice 
fraction but increases with tube diameter. 

On the other hand, different from the aforementioned flow 
patterns, Hirochi et al. [116] observed three types of flow patterns 
of ice slurry with chipped ice particles. The authors believed that 
the behavior of ice slurry was different from that of ordinary 
solid-liquid two-phase slurry due to the cohesive nature of ice 
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particles. The cohesive nature of ice particles dominated ice slurry 
flow pattern while the gravity dominated that without cohesive 
nature. Therefore, the states of moving bed and stationary bed 
were not observed but instead of cluster flow and column flow, 
where the pressure drop also decreased as the reduction of the 
flow velocity, as shown in Fig. 15. 
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Fig. 14. Deposition velocities of ice slurry flow in tubes, (a) Kitanovski [114] and 
(b) Rached et al. [115]. 
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Fig. 13. Ice fraction profiles of ice slurry flow (simulated), base liquid was 10% ethanol/water mixture, d,,=1 mm, D;=27.2 mm, (a) x,=10 wt% and (b) x,=25 wt% [113]. 
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3.2. Rheology of PCS and the friction factor correlations 


3.2.1. Ice slurry 

It was commonly considered that ice slurry with low ice 
fraction, i.e. 6-15% [108], behaved as Newtonian fluid. Therefore 
the flow friction can be easily predicted by empirical correlations 
while the apparent viscosity was usually calculated by Thomas 
relation [108,109] 
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Fig. 15. Flow patterns and pressure drops of ice slurry with chipped ice particles 
in tube [116]. (a) Dispersed flow, U=2.5 mjs, (b) cluster flow, U=1.5 m/s and 
(c) column flow, U=1.0 m/s. 
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Using Poiseuille equation for laminar flow and Blasius correlation 
for turbulent flow, Grozdek et al. [117] concluded that the 
calculation discrepancies were smaller than 10% for both the 
laminar flow of ice slurry with x, < 20 wt% in 9 mm tube and the 
turbulent flow with x,<15wt% in 9mm and 15mm tubes. 
Fernández-Seara et al. [118] also treated the 5-15 wt% ice slurry 
as Newtonian fluid in a fan-coil unit and predicted the flow friction 
with relative error of — 21.5 to 29 % in laminar region by Poiseuille 
equation, while with relative errors of — 6.3 to 9% and — 17 to 1% in 
turbulent region by Petukhov and Blasius equations, respectively. 

However, many other researchers considered ice slurry as non- 
Newtonian fluid and developed semi-empirical correlations based 
on rheological properties, e.g., the Bingham model, the Ostwald 
model (also known as the Power-law model), the Casson model 
and the Herschel-Bulkley model. 

Based on the experimental flow curves of 0-30 wt% ice slurries 
in 10-20 mm tubes as shown in Fig. 16(a), Niezgoda-Żelasko and 
Zalewski [119] determined ice slurry as Bingham fluid, which 
indicated that the ice slurries would behave as Newtonian fluid 
after conquering the yield stress. The corresponding yield stress 
and apparent viscosity for ice slurries with different ice fractions 
were measured, then the friction of laminar flow was estimated by 
Buckingham-Reiner equation, i.e., Eq. (2), with excellent agree- 
ment and the friction factor of turbulent flow was calculated by 
Tomita equation, i.e., Eq. (3), with the uncertainty better than 10% 


1 fr He Het (2) 
Reg 16 6Reż 3f Re8 
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where Regr was calculated by the equation given by Tomita [119]. 
Different from Niezgoda-Zelasko and Zalewski [119], Kumano 
et al. [120] treated ice slurry as pseudo-plastic fluid based on their 
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Fig. 16. Rheological feature determination of ice slurry, (a) shear stress and shear rate of ice slurry flowing through 16 mm tube [119]; (b) flow curve of ice slurry flowing 
through 10.2 mm tube [120] and (c) apparent viscosity of ice slurry measured by rheometer [121]. 
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measured flow curves of 6-24 wt% ice slurries in 4.3-10.2 mm 
tubes, as shown in Fig. 16(b). The authors utilized the Power-law 
model to illustrate the flow characteristics of ice slurry, and 
obtained the corresponding parameters: flow behavior index n 
and fluid consistency coefficient Kk’, then the modified Reynolds 
number, Reyr, was used to describe the flow of power-law ice 
slurries, as given in the following: 


ie Dvp ( 4n o _ 
K \3n+1 D 

It was found that the ice slurry deviated from Newtonian fluid 

more with higher ice fraction and in larger tube. 

Stokes et al. [121] measured the apparent viscosities of 
0-29 vol% ice slurries by a rheometer, results are shown in 
Fig. 16(c). As shown in the figure, there was a large zero-shear 
viscosity at low applied stress, after which the apparent viscosity 
dropped drastically as the increasing shear stress. Therefore, the 
authors suggested that the Herschel-Bulkley model should be 
applied to ice slurries. 

Grozdek et al. [117] tried several types of models and friction 
correlations on the flow of 0-30 wt% ice slurry in 9-25 mm tubes 
and revealed that: Buckingham-Reiner equation for laminar flow 
predicted the experimental results with the relative error within 
+15% for only 60% of all results, and the inclusion of stratification 
effect would extend this error range to 82% of all experimental 
data; Dodge—Metzner correlation for Power-law fluid, i.e., Eq. (5) 
and Tomita equation for Bingham fluid showed good agreement 
with the experimental results in turbulent flow 
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However, Reynolds numbers utilized in the study were not 
corresponding to each rheological model, but just used Reynolds 
number defined as Newtonian fluid (apparent viscosity was 
calculated by Thomas relation). 

Niezgoda-Zelasko and Zelasko [122] proposed a general 
method to estimate the flow friction of ice slurry flowing through 
passages with arbitrary cross-section. The calculation process was 
similar to that of Power-law method, however, the geometrical 
parameters of the flow passage was added to the expressions of 
flow behavior index and fluid consistency coefficient. Reyp was 
then obtained, using which Poiseuille equation for laminar flow 
and Dodge-Metzner correlation as well as Blasius equation for 
turbulent flow could be employed by non-Newtonian ice slurry 
flow. The comparison between the calculated and measured 
friction factors of ice slurry flow through circular tube (16 mm 
in diameter) and rectangular channel (3 mm x 38.5mm and 
78 mm x 26.5 mm) showed the uncertainties were less than 5% 
in laminar region and 10% in turbulent region. 

On the basis of the experimental results obtained by Bédécarrats 
et al. [123] and Grozdek et al. [117], Monteiro and Bansal [124] 
evaluated the mentioned four rheological models. Among these 
models, the Herschel-Bulkley model was considered as the best 
one and the only one which predicted the trend of experimental 
results well due to the consideration of ice particles diameter. 
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3.2.2. MPCS 

It is believed that MPCS can be treated as Newtonian fluid if 
the slurry is homogeneous due to the very fine particles, and this 
conclusion has been verified by many measured viscosities using 
rheometer [9-15,125-128], as given in Table 1, except Zhang 
et al. [129] who obtained pseudo-plastic feature of 10-40 wt% 
MPCS using a rotating cylinder type of rheometer. 

The results reveal the viscosity of MPCS is related with the 
particle fraction and temperature. Vand’s equation was used by 


many researchers [9,10,12,125,127] to predicted the viscosity of 
MPCS, given by the following equation: 


= Up(1—p—Ag?y 7° (6) 


where A is the constant that concerning the particle size, shape 
and rigidity, which is 1.16 with the assumption that the particles 
are ideal rigid sphere. Dammel and Stephan [125] found this 
equation underestimated the viscosity of MPCS, since the parti- 
cles in MPCS was obviously not rigid sphere; Yamagishi et al. [9] 
used A=3.7 while Wang and Lin [127] determined the value 
of 4.4 according to their experimental result. Considering as 
Newtonian fluid, the calculated friction factor using Poiseuille 
equation and Blasius equation given by Yamagishi et al. [9] 
showed perfect agreement with the measured result, as shown 
in Fig. 17(a); the friction factor of MPCS in turbulent flow region 
measured by Wang et al. [12] was higher than the predicted value 
by Blasius equation, as shown in Fig. 17(b), and the authors 
attributed this higher value to the high tube wall roughness; 
Inaba et al. [10] reported the overestimation of Blasius equation 
in MPCS turbulent flow, as shown in Fig. 17(c), the authors 
explained it by the flow drag reduction effect of the fine particles 
(1.5 um) and found that the increase of the ratio of larger particles 
(17 um) caused the smaller discrepancy between the measured 
friction and the calculated values by Blasius equation. 


3.2.3. CHS 

So far, CHS was determined as non-Newtonian fluid in all the 
open literature [24-33,36,37,131]. Based on the measured flow 
curves of TBAB CHS in 16mm tubes, as shown in Fig. 18(a), 
Darbouret et al. [25] believed that TBAB CHS behaved as Bingham 
fluid. The authors presented the corresponding yield stress and 
viscosity which were about 0-14.18 Pa and 0.008-0.17 Pas, 
respectively, depending on the volume fraction of CHS, and found 
good agreement between the measured pressure drop in laminar 
flow region and the predicted values using Buckingham-Reiner 
equation. In addition, Song et al. [27] also treated TBAB CHS as 
Bingham fluid in their study of heat transfer characteristics of the 
slurry. 

Other researchers determined TBAB CHS as pseudo-plastic 
fluid by the flow test in pipes, including Hayashi et al. [24], Xiao 
et al. [26], Ma et al. [28], and Kumano et al. [30]. The results of the 
relation between In(DAP/4L) and In(8v/D) of type B TBAB CHS 
measured by Ma et al. [28] are shown in Fig. 18(b), which 
indicated that the Power-law model was appropriate to TBAB 
CHS. The flow behavior index n and fluid consistency coefficient K’ 
obtained by different authors are summarized and presented in 
Fig. 19, and it was commonly accepted that larger slurry fraction 
caused the slurry deviated Newtonian fluid more. Using the 
modified Reynolds number (Eq. (4)), Ma et al. [28] calculated 
the slurry pressure drop in straight tubes: Poiseuille equation was 
used for laminar flow, and the relative error was + 10% and + 15% 
in 6mm and 14mm tube, respectively; while Dodge-Metzner 
correlation was used for turbulent flow, and the relative error was 
—20% to —25% and +8%, respectively. 

Based on the tested viscosities using rheometer, Suzuki et al. 
[36,37] obtained the shearing-thinning pseudo-plastic behavior of 
TME CHS, and Hashimoto et al. [131] obtained the similar 
behavior of TBAB CHS and TBAF CHS. The latter study also found 
that the TBAB and TBAF solution were both Newtonian fluid when 
the shear rate was low but their apparent viscosity increased 


when shear rate was larger than 500 s~'. 


3.3. Pressure drops and flow friction of PCS 


Compared to single-phase fluid, pressure drops of the two- 
phase PCS are dominated by many other parameters in addition 
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Table 1 
Measured viscosities of MPCS. 
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Authors Materials (including Slurries Temperature, Viscosity, mPas Descriptions 
PCM material and eC 
coating material) 
Yamagishi n-Tetradecane; 3-1000 with 0.2-100 s~! shear 4 decreased with increasing shear rate 
et al. [130] melamine resin . rate 
n-Tetradecane; dp=5 um; A 25 1-7 4 was independent of shear rate 
melamine resin; 1 wt% — ?s=20-30 vol% 
surfactant 
Yamagishi Octadecane; melamine- dp=2-10 um; 25 1.17-10.4 4 increased with particle fraction; nearly 
et al. [9] formaldehyde Ps=7-30 vol% independent of temperature 
Alvarado n-Tetradecane; gelatin dp=2-10 um; 2-12 The ratio to that of water about 4 increased with particle fraction, and was 
et al. [11] X;=4.6-17.7 wt% 1-3.5 independent of temperature 
Wang 1-Bromohexadecane dpa=10.112 um; 20 1.57-8.42 H increased with particle fraction 
et al. amino plastics xs=5-27.6 wt% 
[12,13] 
Inaba et al. Small particle: dpa=1.5 pm; 
[10] n-tetradecane Xs=20 wt%; 
Large particle: dpa=17 um; EF ges increased with particle fraction and decreased 
n-docosane X,=0-50 wt% i G with temperature 
(compared to 
small particles) 
Chen et al. Urea-formaldehyde dp=1-20 um; 10, 20 1.5-4.5 
[14] xs=5-15.8 wt% 
Dammel n-Eicosane; dp=1.5-12 um; 30-60 0.5-4 
and polymethylmethacrylate x,=10-20 wt% 
Stephan 
[125] 
Ho etal. n-Eicosane; urea- dp=4-10 um; 10-50 0.5-2.1 
[126] formaldehyde Xs=2-10 wt% 
Wang and n-Octadecane; d,=0.3-3 um; 20-40 1.5-6 
Lin [127] melamine- xs=5-20 wt% 
formaldehyde 
Zeng et al. 1-Bromohexadecane dp=1-20 um 14.3 10 °C: 5.0 wt% slurry: 4 decreased with the temperature rising 
[128] (C,6H33Br); amino (8.2 um average); 2.05,10.0 wt% slurry: 2.23, 
plastics X,=5.0 wt%, 15.8 wt% slurry: 4.18 
10.0 wt%, 15.8 wt% 
Zhang Nonadecane melamine d,=0-4.5 um; 10-50 10-30 wt% slurry: 2-12; 40 wt% u increased with particle fraction and decreased 
et al. [129] resin X;= 10-40 wt% slurry: 12-35 with temperature; however, the shear stress did not 
vary linearly with shear rate, therefore the fluid was 
considered as pseudo-plastic fluid 
Tumuluri Octadecane gelatin dy=2-10 um; xs=7, 17-32 7 wt% slurry: 1.8-2.2 times as 4 was almost constant with respect to temperature; 
et al. [15] 11 wt% large as that of water; 11 wt% H increased with particle fraction; MPCS was 
slurry: 2.5-3.0 times as large as Newtonian fluid up to 18 wt% fraction 
that of water 
MPCS with multi-walled 11 wt% MPCM; 25-35 10-1000 with 0.2-100 s~' shear u increased with temperature; fluid exhibited non- 


carbon nanotubes 


1.1 wt% MWCNT 


rate 


Newtonian behavior, and Herschel-Buckley model 
was used 


to the flow velocity and tube diameter, including particle fraction, 
particle diameter and shape and so on. Therefore, the former 
measured results of PCS pressure drops exhibit irregularities. The 
recent experimental results of PCS pressure drops are summar- 
ized in Tables 2-4. 


3.3.1. Laminarizing effect 
Most of the experimental results showed that the pressure 
drops of ice slurry were higher than that of the carrier fluid at the 


same flow velocity in tubes [108]. However, Knodel et al. [132] 
found that the flow friction of ice slurry with 2-3 mm ice crystals 
in a 24 mm stainless steel tube reduced with increasing ice fraction 
in a nonlinear manner. The friction factor was close to the value of 
water when the ice fraction was low, and a sudden decline 
occurred as the ice fraction increasing from 2 to 4 wt%, thereafter 
the friction factor changed to another constant value. The authors 
characterized the decline as re-laminarization phenomenon by 
considering the flow turned from turbulent state to laminar state 
as ice fraction increasing. This phenomenon can be explained by 
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Fig. 17. Friction factor prediction of MPCS, (a) Yamagishi et al. [9], 10.1 mm tube; (b) Wang et al. [12], 4 mm tube and (c) Inaba et al. [10], 15 mm tube. 


the laminarizing effect of the ice crystals: the interactions of ice 
crystals and the carrier fluid reduce the fluid turbulence level. By 
comparing the results with their former experiments of ice slurry 
in larger tube diameter (48 and 100 mm tube diameter, 8 wt% 
transition ice fraction), they believed that the interactions between 
ice crystals and tube wall were enhanced when the tube diameter 
was small, and the transitional ice fraction was low. 

The re-laminarization phenomenon was also observed when 
ice slurry flowing through 10-20 mm copper tubes by Niezgoda- 
Żelasko and Zalewski [119] and in 22 mm inner tube of double- 
tube heat exchanger with smooth and corrugated surfaces by 
Bédécarrats et al. [123]. As shown in Fig. 20, ice slurry flowing at a 
certain velocity is gradually damped from fully turbulent to 
laminar state as the ice fraction increasing. In the fully turbulent 
region, the pressure drop is nearly constant when the ice fraction 
is low, and then gradually increases with ice fraction. Bédécarrats 
et al. [123] believed that the liquid phase dominated the flow 
pressure drop when the ice fraction was low, therefore the 
pressure drop did not change significantly in the turbulent region. 
The pressure drop increases more rapidly in the laminar region 
than that in turbulent region, which was attributed by Bédécar- 
rats et al. [123] to the increasing friction of ice particles on the 
tube wall. Moreover, pressure drop reduces with ice fraction 
when ice slurry flows in the transitional region, and the transi- 
tional ice fraction moves towards the high ice fraction with 
increasing flow velocity as clearly shown in the Fig. 20(a). 

Similar phenomena were observed in the flows of MPCS and 
TBAB CHS. Yamagishi et al. [9] noticed the lower pressure drop of 


30 vol% MPCS at the flow velocity of 2.0-2.5 m s7! compared to 
that of water, and attributed the laminarizing effect to the 
depression of longitudinal eddy diffusivity of slurry since the 
particle size was smaller than the turbulence length. The authors 
also pointed the laminarization was beneficial for pumping power 
reduction. Xiao et al. [142] and Ma et al. [28] reported the 
laminarization during the flowing of TBAB CHS in tubes, the 
variations of pressure drop with particle fraction at the certain flow 
velocity were quite similar to that of ice slurry, as shown in 
Fig. 21(a) and (b), the reduction of pressure drop were both noticed 
occurred at the slurry fraction range of 12-16 vol%. Ma et al. [28] 
also noticed that the laminarizing effect was particle fraction and 
flow velocity dependent, which indicated the laminarizing effect 
would disappear when the particle fraction and flow velocity were 
both high. This was caused by the enlarged particle interactions and 
therefore the mechanical friction was larger. 


3.3.2. Influence of flow passage geometry on flow friction 

Different from that of the single phase fluid, the flow friction of 
solid-liquid two-phase slurry is composed of viscous friction 
caused by the fluid viscosity and the mechanical friction which 
is resulted from the interactions between solid particle and the 
flow boundary [144]. Therefore, the flow friction of PCS is largely 
related to both the slurry fraction and flow passage geometry, 
such as passage diameter, roughness and pattern. 

Relationships between the flow friction and Reynolds number 
of ice slurry in tubes with different diameters were obtained by 
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Gozdek et al. [117], Niezgoda-Zelasko and Zalewski [119] and 
Illán and Viedma [133]. As shown in Fig. 22(a)-(c), the experi- 
mental results obtained by these studies show similar friction 
factor-Reynolds numbers relationships though the used Reynolds 
numbers were defined in different ways: Grozdek et al. [117] 
determined the Reynolds number using the viscosity calculated 
by Thomas relation (Re;,), Niezgoda-Zelasko and Zalewski [119] 
used the Reynolds number based on the Bingham model (Rep), 
while Illán and Viedma [133] employed the referenced Reynolds 
number (Reqs) corresponding to the carrier fluid at 0 °C. In the 
laminar flow region, it was noticed that both the higher ice 
fraction and the larger tube diameter resulted in larger flow 
friction, and the corresponding deviation from that of the 
Newtonian fluid also became larger. Grozdek et al. [117] illumi- 
nated the influence of the tube diameter on the ice slurry flow: ice 
slurry flowing through smaller tubes tended to be more homo- 
geneous at larger mean velocity but still in laminar state, while 
the mean velocity was relatively low in larger tube which resulted 
in heterogeneous flow with more ice particles at the tube top but 
less at the bottom, therefore the distorted flow pattern led to 
higher flow friction. Kumano et al. [120] observed the similar 
influence of tube diameter on the flow of ice slurry while narrow 
tubes with the tube diameter of 4.3 mm, 7.5 mm and 10.2 mm. 
Ma et al. [28] regarded that the particle-tube wall interac- 
tions was enlarged when TBAB CHS flowing through tube with 
smaller diameter, hence the fluidity of the slurry deviated from 
Newtonian fluid much more compared to that of the slurry 
flowing in larger tube, and then the flow friction was higher 
due to the larger mechanical friction loss in the slurry. In their 
study, the measured flow friction factor in 6 mm tube was about 
20-25% larger than the predicted value by Dodge-—Metzner 
correlation as mentioned before, while the prediction error in 
the case of 14mm tube was only +8% larger. They believed this 
tube wall constraint effect on slurry flow would disappear if the 
tube diameter was large enough compared to the particle size. 
Both Bédécarrats et al. [123] and Illán and Viedma [133] 
observed reasonably higher pressure drop of ice slurry in corru- 
gated tube compared to that in smooth tube. In the former study, 
different from the slight increase in smooth tube, the pressure 
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fraction in [26] is volume fraction while in [28,31] are mass fraction). 


Table 2 


Pressure drops of ice slurries. 


Authors Base liquids Slurries Geometries Flow rates Correlations Descriptions 
Knodel et al. Water d,=2-3 mm; Stainless steel tube, 2.8-5.0 m s7! f =0.946 x 0.184Re~°?, 38,000 < Re < 74,000 AP/L was nearly constant when x, < 2 wt% and x, > 5 wt%; 
[132] Xs < 12 wt% Dj=24 mm, L=4.596 m AP/L reduced when 2 <x, <5 wt% due to the laminarizing 
effect 
Meewisse 5-39% Various d,=0.1-1 mm; Circular tube, D;=20 mm 3ms-! AP/L increased with ice fraction; AP/L depended on the 
and infante water-brine Xs < 30 wt% added freezing point depressant. Slurry with methanol, 
Ferreira NaCl and potassium formate had lowest pressure drop 
[134] 
Bellas et al. 5% Propylene  d,=0.17- Plate heat exchanger, 24 1.0-3.7 m? ho! AP increased with v, and xs; only a 3-8% increment in AP 
[135] solution 0.6 mm; plates, Dı=4 mm, was observed as x, varying from 5.0 to 20.0 wt% 
Xs < 20 wt% W=112 mm, H=310 mm 
Nørgaard 16% Propylene x; <30 wt% Plate heat exchanger, 10 plates, 50-300 L h~! AP of ice slurry in heat exchange generally increased with v, 
et al. [136] glycol solution Dn=4 mm, W=70 mm, and xs; slight reduction was observed when flow rate 
H=300 mm increased from 100 to 150 Lh~! with inlet 30 wt% ice 
fraction 
Hägg [137] 33.2 wt% Xs < 30 wt% Stainless steel tube, D;=9, 15, 0.5-2 m s7! AP/L increased with vs and xs and with reduced pipe 
HYCOOL 25 mm, L=1.2 m dimensions 
solution and 
35 wt% ethanol 
solution 
Niezgoda- 10.6% Ethyl dp=0.05- Copper tubes, D;=10, 16, 0.1-4.5ms7! Buckingham-Reiner equation for laminar flow, excellent APJL increased at the fully turbulent and more rapidly at the 
Zelasko and alcohol 0.35 mm, 20 mm, L=4.6 m agreement; Tomita method for turbulent flow, uncertainty laminar flow regions with x, at a certain vs, but decreased at 
Zalewski solution Xs < 30 wt% was better than 10% the transitional region, while the Xtra increased with vs; Retra 
[119] increased with xs 
Shire et al. 4.76 wt% NaCl 0.02-0.05 mm in Borosilicate glassware, 0.5-2 m s7! AP/L increased with v, and x,; AP/L depended largely on a 
[138] brine width, 0.1- D;=25 mm, L=1 m critical yield stress for high xs 
0.2 mm in length, 
41 <x, < 67 wt% 
Illán and 9% NaCl brine d, {~ 0.5 mm, Stainless steel corrugated tube, 700, f= 6.12K°*6Re P16, x is the severity number of the AP was almost constant with x, 
Viedma Xs < 25 wt% D;=16 mm, L=1.25 m 1200 kg h“! corrugated tube 
[139] 
Shire etal. 5% NaCl brine 0.02-0.05 mm in Plate heat exchanger, 12 0-2.5 Ls! AP increased with v, and x,; AP had sharp increase when x, 
[140] width, 0.1- plates, D,=7 mm, W=250 m, was higher than 40 vol% in the tubular heat exchanger; 
0.2 mm in length, H=1000 m, tubular heat required pressure drop was lower than the predicted from 
10 <xs < 60 vol% exchanger, D,=24 mm, L=3 m models; mixing mitigated the block caused by high ice 
fraction as well as large ice particles 
Bédécarrats 10% Ethanol dp x 0.1 mm, Smooth and corrugated tube, 0.3-1.9 m s7! At a certain vs, AP/L increased slowly with x, within smooth 
et al. [123] solution Xs < 30 wt% Dj=22 mm, L=2 m tube, but the increment stopped at critical x,, while AP/L 
was nearly constant within corrugated tube; then AP/L 
reduced to a minimum value and then rapid increase; Xtra 
increased with v, 
Grozdek 10.3 wt% dp x 0.3 mm, Stainless steel tube, D;=9, 15, 0-2.6 m s7! The authors tried several kinds of correlations, the best ones Ice slurry flow behaviour was time dependant; AP/L generally 
etal. [117] Alcohol Xs < 30 wt% 25 mm, L=1.2 m were the modified Buckingham-Reiner for laminar flow increased with x, at a certain v, however, the delayed 
solution and Dodge-Metzner and Tomita for turbulent flow laminar-turbulent transition with higher x, resulted in lower 
APJL than that of carrier fluid sometimes; when x, < 10- 
15 wt%, the slurry could be treated as Newtonian fluid 
9% NaCl brine d,~0.5 mm, Stainless steel tube, 0.2-5 m s7! Smooth tube (subscript cf stands for the carrier liquidat f increased with x, in laminar region, but was nearly 
Xs < 26 wt% dm/Di=0.014-0.031 0 °C): Laminar flow: f = ARe‘% constant in turbulent flow region; f increased with the 
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Table 2 (continued ) 


Authors Base liquids Slurries Geometries Flow rates Correlations Descriptions 
Illán and A=0.00024(1—0.1x, /100)B!*!”2 decrease of dm/D; in laminar region, but was nearly constant 
Viedma 4,\ 0.578] 4 in turbulent flow region; sharp increase was found when xs 
[133] x [2:95(t) | exp((a+1)b) was up to 26 wt% in corrugated tube 
a = (0.1427—0.174 In B), 
b = (0.22435—0.000308p,,)xs/100 
1.166 
B=1 .295D? pe” 100% usp (8) e 
Turbulent flow: f = 0.316cRe;°?94, 
c = 1.4074(1—0.1x; /100)e— 9.97 185%s 
d=1.1—0.02847x, 
Kumano 5wt% Ethanol d,=0.1-04mm, Dj=4.3, 7.5, 10.2 mm, L=2.5 m Homogeneous flow without ice particles accumulated at the 
et al. [120] solution Xs < 25 wt% tube top; AP/L increased with x, in the laminar flow region; 
APJL was constant or decreased slightly with x, in the turbulent 
flow region when x; was low, but increased with xs when xs 
was high; slurry flowing through larger tube deviated more 
from the Newtonian fluid and had higher flow friction 
Fernández- 10 wt% dp=0.025- Fan coil heat exchanger, 0.3-1.1 ms! Poiseuille model for laminar flow with relative errors of- AP/L increased with v, and xs 
Seara etal. Ethylene glycol 0.25 mm, smooth copper tube, Dj=3/8” 21.5% to 29%, Petukhov equation for turbulent flow with 
[118] solution Xs < 20 wt% relative errors of -6.3% to 9%, Blasius equation for turbulent 
flow with relative errors of -17% to 1% 
Table 3 
Pressure drops of MPCS. 
Authors PCM and shell material Slurries Geometries Flow rates Correlations Descriptions 
Yamagishi Octadecane; melamine- d,=2-10 um, Stainless steel pipe, 0-3.0 m s7! Poiseuille equation for laminar APJL increased with vs; AP/L of 7-20 vol% slurry was higher than that of water 
et al. [9] formaldehyde @s=7-30 vol% D;=10.1 mm, L=800 m flow and Blasius equation for and increased with @,; laminarization was observed when ps=30 vol%, which 
turbulent flow with perfect caused APL smaller than that of water at v,=2.0-2.5 m s7! 
agreement 
Rao et al. n-Octadecane; d,=1-5 um, Rectangular copper 0.1-1.4 m s7! AP increased with v, and xs 
[141] polymethylmethacrylate x,=5-20 wt% minichannels, 
2 mm x 4.2 mm x 150 mm 
Alvarado n-Tetradecane; gelatin dp=2-10 um, Copper pipe, Dj=10.9 0.5-1.5 m s7! in AP/L increased with v,; In 10.9 mm smooth tube, AP/L was slightly higher 
et al. [11] X,=5.9-15.2 wt% (smooth), 8 (smooth), 10.9 mm tube; than that of water (calculated value); In 8 mm tube, AP/L of 5.9 wt% slurry 
8 mm (enhanced), 1.2-2.5ms~! in was smaller than that of water (calculated value) due to the turbulence 
L=12.2m 8 mm tube suppress effect by small particles; AP/L in smooth and enhanced tubes were 
similar 
Chen et al. 1-Bromohexadecane; dp=1-20 um, Stainless steel pipe, 0.2-1.0 m s7! AP increased with v, and increased slightly with temperature; AP was 
[14] urea-formaldehyde xs=15.8 wt% at D;=4 mm, L=1.85 m higher than that of water at the same vs, but the difference reduced with 
10 and 20 °C increasing flow velocity 
Dammel n-Eicosane; d,=1.5-12 um; Rectangular copper 0.2-0.8 kg min-! AP increased with v, and x, 
and polymethylmethacrylate x,=10-20 wt% minichannels, 
Stephan 2mm x 4.2 mm x 390 mm 
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Pressure drops of CHS. 
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Authors Initial Slurries Geometries Flow rates Correlations Descriptions 
solution 
concentration 
Darbouret 20 wt%, TBAB CHS, Type Dj=16 mm, L=1 m 0-300 Lh7! Buckingham-Reiner AP increased with xs, and increased linearly 
etal. [25] 35 wt% A: Xs=4- equation for laminar with flow rate; Bingham fluid was 
55 wt%; Type B: flow determined for TBAB CHS 
X,;=28-43 wt% 
Xiao et al. 22 wt%, TBAB CHS: PVC pipe, D;=21 mm, L=2m 0-5.25 m s7! Laminar flow: AP increased with vs; AP sometimes 
[26] 30 wt% @s=0- Poiseuille equation; decreased with ps; power-law feature was 
16.24 vol% Turbulent flow: determined for TBAB CHS 
fr=0.3226Remn*?”! 
Suzuki 25 wt% TME CHS: Dj=5 mm, L=0.5 m Re (calculated Drag reducing by the used surfactants was 
et al. X,=0-18.8 wt%; based on water effective for TME CHS when the flow rate and 
[36,37] Ethoquad 0/12 viscosity) was in the slurry fraction were both high 
as surfactant the range of 
2000-20,000 
Xiao et al. Type A TBAB PVC pipe, D;=16 mm, L=2.8m 0-3.5 m s~! AP increased with v,; as the increase of x,, 
[142] CHS: x,=0- the flow of CHS was in the turbulent region 
27 wt% first and AP increased slightly; when x, was 
in the range of 10-15 wt%, laminarization 
was observed, and AP dropped largely; then 
AP increased rapidly in the laminar flow 
region 
Ma et al. 10 wt% Type B TBAB Copper tube, Dj=6, 14 mm, 0-2.5ms~'in Laminar flow: In laminar region, AP increased slightly with 
[28] CHS: g,=0- L=1.5m 6mm pipe, 0- Poiseuille equation; ọsin 6 mm tube but changed little in 14 mm 
20 vol% 1.2ms™! in Turbulent flow: tube; in the transitional and turbulent 
14 mm tube Dodge-Metzner regions, the laminarizing effect caused a 
correlation decrease of AP with ps, however, the 
laminarizing effect was volume fraction and 
flow velocity dependent 
Kumano 22.5 wt% Type A TBAB AP increased with v,; AP increased with x, in 
et al. [30] CHS: x;=0- laminar flow region but slightly decreased in 
20 wt% Stainless steel pipe, Dj=4.3, turbulent flow region; In turbulent flow 
8.5 wt% Type B TBAB 7.5, 10.2 mm, L=1 m region, the flow characteristics changed at a 
CHS: xs=0- certain flow velocity while xs was about 
20 wt% 15 wt%, then AP dropped a lot 
Hashimoto 20 wt% TBAB TBAB CHS: The addition of surfactant caused the 
et al. [131] solution; xs=0-70 wt%; drastically reduction of the viscosities of both 
18 wt% TBAF TBAF CHS: TBAB CHS and TBAF CHS 
solution Xs=0-42 wt%; 
Ethoquad 0/12 
as surfactant 
Ma et al. 10 wt% Type B TBAB Plate heat exchanger, 10 2.5-13 Lmin“! f=51.002Remg”, AP increased with vs; AP was always higher 
[33] CHS: ps=0- plates, D»=4 mm, W=72 mm, 60 < Remr < 400 than water; AP was almost constant when 
17.5 vol% H=310 mm 5< P; < 13.6 vol%, but increased a lot when 
@s=17.5 vol% 
Ma et al. 15 wt% Type B TBAB Copper pipe fittings, 14 mm 1-8ms—lin The local loss friction increased slightly with 
[29] CHS: x,=0- 90° elbow, 14-6 mm sudden 6mm pipe, x, when x, was low but increased rapidly if x, 
30 wt% contraction, 6-14 mm sudden 0.4-2.4 m s™! in was larger than 20 wt% 
expansion, 14 mm tee 14 mm tube 
Suziki et al. 35 wt% Ammonium Copper tube, D;=13 mm, 0.5-1.5 m s7! AP of slurry was slightly larger than that of 
[143] ammonium alum hydrate L=2 m water at the same Re. surfactant was 
alum slurry: effective, with which AP was even lower 
Xs=12 wt% than that of water 


drop of ice slurry flowing in the turbulent region in corrugated 
tube was nearly constant or only had slight reduction as increas- 
ing the ice fraction. In the latter study, pressure drop in the 
laminar region was noticed a sharp increase when the ice fraction 
was up to 26 wt%, and the authors also successfully introduced 
the severity parameter relevant to the corrugated geometrical 
characteristics to predict the pressure drop. Alvarado et al. [11] 
compared the pressure drops of 5.9 wt% MPCS in regular and 
enhanced 8 mm tubes (the tube roughness was considered as 
0.1 mm), however, the results showed the similar pressure drop 
in magnitude. 


3.4. Heat transfer characteristics of PCS 


Numerous experiments have been performed to investigate 
the overall and local heat transfer coefficients of PCS flowing 
through the heated tube and heat exchangers, the relating investi- 
gations were summarized in Tables 5-7. It was expected that the 
latent heat of the solid particles in slurry could improve the heat 
transfer of ice slurry, and many reported investigations validated 
this improvement, as shown in Fig. 23 which presented the heat 
transfer enhancement of ice slurry summarized by Stamatiou and 
Kawaji [112], however, the comparison may conducted under the 
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Fig. 20. Pressure drops of ice slurry as the increase of ice fraction, (a) 16mm 
copper tubes [119] and (b) 22 mm inner tube of double-tube heat exchanger with 
smooth surface [123]. 


same Reynolds numbers since the flow turbulence of slurry was 
damped by the particle flow. Some researchers also attributed 
this heat transfer coefficient improvement to the turbulence in 
the near-wall thermal boundary layer due to the movement of 
solid particles. However, it was believed that particle fraction, 
flow velocity, heat flux, particle shape and size, pipe shape and 
dimension, heat exchanger type and size imposed influences on 
the heat transfer of ice slurry, therefore the combined effect can 
be complicated. 


3.4.1. Local heat transfer along the flow direction 

Stamatiou and Kawaji [112] obtained the variations of local 
Nusselt numbers of ice slurry flowing though heated rectangular 
channels, as shown in Fig. 24(a). Nujocq, Was large at the channel 
inlet and decreased slightly towards the exit of the channel. The 
authors attributed the large Nujo-q; at the inlet to the developing 
flow and the following decrease to the more developed flow. Hagg 
[137] reported the similar slight reduction of local heat transfer 
coefficient of ice slurry along the flow direction in 21mm 
stainless steel pipe, as shown in Fig. 24(b). The authors believed 
that the phase change and turbulence caused by ice particles in 
the thermal boundary layer led to smaller temperature gradient 
at the entrance and therefore high hican and the turbulence 
reduced as the ice particles melting along the flow direction and 
therefore hiocaı decreased. Moreover, Ionescu et al. [147] reported 
hiocaı Of ice slurry decreased greatly at the channel inlet, followed 
by nearly constant values. 
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Fig. 21. Pressure drops of TBAB CHS as the increase of particle fraction, (a) type A 
TBAB CHS in 27 mm PVC tube [142] and (b) type B TBAB CHS in 6 mm copper tube 
[28]. 


Yamagishi et al. [9], Alvarado et al. [11], Wang et al. [12,13] 
and Tumuluri et al. [15] obtained similar hj.,q; variations of 
turbulent flow of MPCS in heated tubes, typical results were 
shown in Fig. 25. The local heat transfer coefficient increased at 
the tube inlet, where the slurry temperature was lower than its 
phase change temperature (supercooled slurry). The increasing 
Niocat Was Caused by the gradually increasing melted particles near 
the tube wall; then the slurry temperature reached its phase 
change temperature, and hoca decreased in this region since the 
number of solid particles decreased gradually and the fully 
melted layer near the tube wall deteriorated the heat transfer 
from tube wall to fluid core. After the slurry was finally fully 
melted, Ajocai increased again due to the rising slurry temperature. 
Yamagishi et al. [9] and Wang et al. [12,13] also showed the 
variation of hoca in laminar flow region, which generally 
decreased along the flow direction. The former study pointed 
that the phase change particles in the laminar flow was possibly 
not completely melted, leading to generally decreasing solid 
particle number and then decreasing heat transfer coefficient. 

As shown in Fig. 26(a), Kumano et al. [31] reported the 
generally decreasing hiocaı of TBAB CHS in laminar flow region, 
and attributed the decrease to the undeveloped boundary layer of 
the slurry. The authors also found that the variation of Njocq along 
the heating tube changed little under the turbulent flow condi- 
tion, which was caused by the rapidly developed boundary layer. 
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Fig. 22. Relationships between Reynolds number and friction factor of ice slurries, (a1) and (a2) Niezgoda-Zelasko and Zalewski [119] and (b) Grozdek et al. [117]; (c1) and 


(c2) Illán and Viedma [133]. 


The results reported by Ma et al. [28] were shown in Fig. 26(b), 
and it can be seen from the figure that hiocaı decreased along the 
heating tube but increased after all the solid particles were fully 
melted. This variation was explained by the similar reason to that 
of MPCS, except there was no supercooling region for TBAB CHS. 


3.4.2. Influence of particle fraction on heat transfer 

Heat transfer performances of ice slurry flowing through 
heated circular tubes have been investigated by many research- 
ers. However, sometimes quite different even opposite character- 
istics of the heat transfer coefficients with ice fraction were 
reported, and the same phenomena happened to MPCS and CHS. 

Similar to the re-laminarization phenomenon of the ice slurry 
flow as aforementioned, Knodel et al. [132] also observed the 
reduction of the Nusselt number of ice slurry in the heated tube, 
as shown in Fig. 27(a). This reduction was attributed to the 


damped flow turbulent degree by laminarzing effect. As the 
reported by Niezgoda-Żelasko [145] and Niezgoda-Żelasko and 
Żelasko [148], there was a slight increment of the heat transfer 
coefficients of ice slurry as ice fraction increasing from 5 to 
30 wt% in the laminar flow region, while the heat transfer 
coefficients were nearly constant in the turbulent flow region, 
as shown in Fig. 27(b). Clearly, there is a sharply decline of the 
heat transfer coefficients as ice slurry flowing in the transitional 
region. Moreover, the authors pointed that the influence of ice 
fraction on heat transfer was more distinct in the laminar flow 
region while the heat transfer coefficients were more strongly 
dominated by flow velocity than by ice fraction in the turbulent 
flow region. The same conclusion was also obtained by Grozdek 
et al. [149] and Kumano et al. [150] in the case of heated tubes, 
and Bédécarrats et al. [123] in the case of double-tube heat 
exchanger. Grozdek et al. [149] observed generally increased heat 
transfer coefficient with ice fraction as well as flow rate, as shown 


Table 5 


Heat transfer of ice slurries. 


Authors 


Knodel et al. 
[132] 


Bellas et al. 
[135] 


Norgaard et al. 
[136] 


Stamatiou and 
Kawaji [112] 


Niezgoda- 
Zelasko [145] 


Lee et al. [146] 


Ionescu et al. 
[147] 


Base liquids 


Water 


5% Propylene solution 


16% Propylene glycol 
solution 


6.2 wt% NaCl brine 


10.6% Ethanol solution 


6.5% Ethylene glycol 
solution 


Slurries 


d,=2-3 mm, 
Xs <12 wt% 


d,=0.17-0.6 mm, 
Xs < 22 wt% 


Xs < 30 wt% 


dp=0.1-0.2 mm, 
Xs <25 wt% 


d,=0.1-0.15 mm, 
Xs < 30 wt% 


dp x 0.27 mm, 
Xs < 26 wt% 


Composed by ice 
particles and oil, 
dp x 0.1 mm, 
Xs=20, 25, 30 wt% 


Geometries 


Stainless steel tube, D;=24 mm, 
L=4.596 m 


Plate heat exchange, 24 plates, 
D =4 mm, W=112 mm, H=310 mm 


Plate heat exchange, 10 plates, 
D =4 mm, W=70 mm, H=300 mm 


Rectangular channels, 
533 mm x 305 mm x 25 mm, 
533 mm x 100 mm x 13 mm 


Copper tubes, D;=10, 16, 20 mm, 
L=4.6 m 


Double-tube heat exchanger, inner 
tube: D;=13.84 mm, outer tube: 
D=25 mm, L=1.5m 


Rectangular channel, 
8 mm x 80 mm x 1000 mm 


Flow rates 


2.8-5.0 m s7! 


1.0-3.7 m? h`! 


50-300 L h7! 


0.1-4.5 m s~! 


800- 
3500 kg m~? s7! 


Correlations 


Num = 0.885NUpetukhov 


Laminar flow: 
Nu) = 4.0Gzef 86x93 (115 / py)?” 


Turbulent flow: Num = Nucn[1 +1.85x 
108x272 Reg (ig / Hw) ] 


Nucn was the single-phase Nusselt 
number calculated by Gnielinski 
equation 


Laminar flow: Nu = 0.78(Re Pr 5)°7 


AH AX —0.027 dp —0.29 
CAT 100 D; 


Turbulent flow: Nu = 0.031(Re Pr)??? 
AH Ax\ 0°18 7,021 
(Garioo) (B) 


ô was the Bingham term 


Num(x) indicated Nu(x) integrated 
from 0 to x 


Laminar flow: 
-for a constant wall temperature: 
—1.12 
Num(x) =0.34 (zp)  Re042 pro33 
a +X5)73 


-for a constant heat flux: 


—1.49 
Num(X) =0.22 (zp)  Re051 pro33 
(d +x)?7 


Transition flow: 


Num(x) = 0.03 (=n) 
(1+xs) 


-0:35 00.75 Pr0-33 


2.45 


Comments 


Num reduced when x, < 6-8 wt% 


Aoverau remained fairly constant as the increase of xs; Noveraut 
increased with vs; 30% increment in heat transfer capacity 
was observed as the x, varying from 5.0 to 20.0 wt% 


Aoveray increased with v, and xs; increment with x, became 
little at high vs 


Nujocat Was high at the inlet and decreased slowly along the 
flow direction; Nuca and Num increased with x, and heat 
flux 


h changed little with heat flux density; h increased slightly 
with x, in laminar region; h decreased with x, at the 
transitional region; h was more strongly dependent on vs 
than on xs in turbulent region; h of 10-20 wt% ice slurry 
can be lower than that of single-phase carrier fluid; Nu was 
much larger than that of carrier fluid 


At low flow rate, Moveran did not change with x, when xs was 
smaller than about 10%, then increased sharply with x,; 
Aoverau increased gradually with x, at high flow rate range; 
hot water temperature had little impact on hoverau 


Aiocat had a great decrease at the inlet, followed by a nearly 
constant value; hjocq increased with v, and x, 


bros 
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Niezgoda- 
Zelasko and 
Zelasko [148] 


Gozdek et al. 
[149] 


Bédécarrats 
et al. [123] 


Illán and 
Viedma [158] 


Illán and 
Viedma [133] 


Fernandez- 
Seara et al. 
[118] 


Kumano et al. 


[150] 


10.6% Ethanol solution 


10.3% Alcohol solution 


10% Ethanol solution 


9% NaCl brine 


9% NaCl brine 


10 wt% Ethylene glycol 
solution 


5 wt% Ethanol solution 


d,=0.1-0.15 mm, 
Xs < 30 wt% 


dp x 0.3 mm, 
Xs < 22 wt% 


dp x 0.1 mm, 
Xs < 30 wt% 


dp x 0.5 mm, 
Xs < 26 wt% 


dp ~ 0.5 mm, 
Xs < 26 wt% 


d,=0.025- 
0.25 mm, 
Xs < 20 wt% 


dp x 0.15 mm, 
Xs <25 wt% 


Copper tubes, D;=10, 16, 20 mm, 
L=4.6 m, rectangular tube, 
0.0078 m x 0.0265 m x 3 m, 
rectangular (slit) tube, 

0.003 m x 0.0358 mx 2m 


Stainless steel tube, D;=21 mm, 
L=0.94m 


Double-tube heat exchanger with 
smooth and corrugated tubes, inner 
tube: D;=22 mm, L=2 m, outer 
tube: Dj=33 mm, L=1.5m 


Fan-coil heat exchanger, cooper tube, 
D;=8.8 mm 


Stainless steel tube, d,,/D;=0.014- 
0.031 


Fan coil heat exchanger, smooth 
copper tube, Dj=3/8” 


D,;=4.3, 7.5, 10.2 mm, L=2.5 m 


0.1-4.5 m s~! 


0-1.25 m s7! 


0.3-1.9 m s7! 


1000-4600 kg 
h-! 


0.2-5 m s7! 


400, 600, 
800 kg h`! 


Turbulent flow: 
-for a constant wall temperature: 


Num(X) = (1.57 +26.1x975) 
Gz?33 + 1.76x2 —1.1x; 


-for a constant heat flux: 


Num(X) = (2.06 + 78.8x! 85) 
Gz?33 + 0.55x2 —0.54x, 


Laminar flow: 


n Pe y 
= m(_AH_ Ax \" (dp\P (Tp 
Nu =A Gz (eir) (3) (eae 


Turbulent flow: 
Nu =A(Re Pr)™ ()'an parameters 
are seen in [148] 


Laminar flow: 
Nu(x) = 87.3G2!(x5/100)°°7 


Turbulent flow: 
My = 1+0.00103x,—2.003 


—0.19230=%s/100 O.339RE pF 


Rey 3 (Xs/100) 10 


Smooth tube: 

Nu = (0.845—0.01216x;)Nug, 

Nug = 0.0344(Re-y—1500)°7*Pr37 
Corrugated tube (x was the corrugated 
parameter): 

Nu = (0.892—0.01824x; )Nucy, 

Nug = 0.37415 (Rep —1 500)°”4Ppro4 


Laminar flow: 


0.822 
Nu =50.1Re9 074 (%5) (2) 


Turbulent flow: Nu = 0.023Rey°8Pr!/3 


Similar to Niezgoda-Zelasko [145]; smaller hydraulic 
diameter of channel caused smaller increase in h of ice 
slurry compared with h of carrier fluid 


h changed little with heat flux; h generally 

increased with x, and v,; x, had greater impact on h than 
vs; h of 10-15 wt% ice slurry had little increase or no 
increase in comparison with carrier fluid; heterogeneous 
flow pattern of ice slurry in laminar region caused higher h 


At a certain vs and as increasing Xs, Rovera increased up to a 
critical x, within smooth tube, followed by a reduction to a 
minimum value and then a rapid increase; Noverat 
decreased first within corrugated tube, then increased; 
Avera depended on heat flux and flow pattern; Rovera Was 
more strongly dependent on vs than x, in turbulent region 


Heat transfer rate increased with vs; frost formation, 
partial blockage and excessive diameter ratio caused the 
worse behaviour prediction 


Num decreased with x, 


Heat transfer rate of ice slurry was higher than chilled 
water by the factor of 3.7-4.6 depending on ice fraction 
and fan rotation rate; 80-88% of the overall thermal 
resistance was on the air side 


h increased with x, and vs; h depended more strongly on x, 
than on vs in laminar flow region; in turbulent region, h 
was nearly constant at low xs but increased largely with vs; 
hħiocaı Changed little with distance 
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Table 6 


Heat transfer of MPCS. 


Authors 


Goel et al. 
[151] 


Yamagishi 
et al. [9] 


Inaba et al. 
[10] 


Alvarado 
et al. [11] 


Wang 
et al. 
[12,13] 


Rao et al. 
[152] 


Diaconu 
et al. 
[153,154] 


Wang and 
Lin [127] 


Wang 
et al. [155] 


Materials 


n-Eicosane, urea- 
formaldehyde 


Octadecane; melamine- 
formaldehyde 


Small particle: n- 
Tetradecane; large 
particle: n-Docosane 


n-Tetradecane; gelatin 


1-Bromohexadecane, 
amino plastics 


n-Octadecane, 
polymethylmethacrylate 


Rubitherm RT6, 
polycyclic cell 


n-Octadecane, 
melamine- 
formaldehyde 


Slurries 


d,=100, 250 um, 
@s=5-20 vol% 


dp=2-10 um, 
@s=7-30 vol% 


dpa=1.5 um, 
Xs=20 wt%; 
dp=17 um, xs=0- 


50 wt% (compared to 


small particles) 


dp=2-10 um, 
Xxs=5.9-15.2 wt% 


dpa=10.112 um, 
Xs=5-27.6 wt% 


dpa=4.97 um, 
Xxs=5-20 wt%; 


dpa=2.24 pm 


d,=0.3-3 um, 
Xxs=5-20 wt% 


dp=0.5-5 um, 
X,=5-20 wt%, with 
0.5 wt% TiO, 
nanoparticles 


Geometries 


Copper tube, D;=3.14 mm, 
L=04m 


Stainless steel tube, D;=10.1 mm, 
L=800 D; 


Stainless steel tube, D;=15 mm, 
L=5.85m 


Copper tube, Dj=10.9 (smooth), 8 
(smooth), 8 mm 
(enhanced)L=12.2 m 


Stainless steel tube, D;=4 mm, 
L=1.46m 


Rectangular copper minichannels, 
2mm x 4.2 mm x 150 mm 


Helically coiled copper tube in 
tank, D;=10 mm, 140 mm helix 
diameter; Water flowed inside tube 
while MPCS was contained in tank 


Aluminum rectangular tube, 
12.8 mm x 1.8 mm x 900.2 mm 


Copper tube, Dj=3.97 mm, 
L=1.834m 


Flow rates 


Reynolds number 
was 200, 1000 
(only for 5 vol% 
slurry) 


0-3.1 m s~! 


300 < Rey < 8000 


400 < Rem < 3500 


0.05- 


0.35 kg min7! 


4.17-10.0 g s~! 


Correlations 


Laminar flow: Nujocat = CNUshan and London» 
c was 1.336-1.418 in the case of 5-15.8 wt% 


MPCS, + 15% relative error 
Num = 0.8148Re2493 Pr9.4836 


Ste~°1277((Ly +. Ly)/D]?3°? Ly +L was the 
phase change region, + 10% relative error 


for 5-27.6 wt% MPCS 
Turbulent flow: Num = 4.8527 x 
1074 Re0.7733 pr2.7941 ste0.3159 

m m 


[da + L2)/D]-9333 (Lp / hy)? 
+ 10% relative error for 5-10 wt% MPCS 


Nature convection correlation, 


Num = 0.802Ra°?78, 10"! < Ra < 1013 


Comments 


Stefan number was the most dominant factor on heat 
transfer; lower Stefan number (or heat flux) caused 
lower wall temperature; effect of ps for a given Stefan 
number could be expected when s was 5 vol%; larger 
particle led to lower wall temperature; the effect of 
homogeneity could be negligible 


In laminar flow region, h decreased along the flow 
direction; in turbulent flow region, h increased in the 
supercooling region along the flow direction, and then 
decreased in the melting region, then increased after 
finishing the melt; higher heating power led to lower h; 
h generally increased with x, and v, 


Nujocat decreased along the flow direction and reached 
nearly constant when fully melted; Num increased with 
vs and the ratio of large particle 


The variation of h along the flow direction was similar to 
the result in [6]; h increased significantly with v,; 
increase of x, caused lower flow turbulence and then 
lower h, and h was lower than that of water at the same 
flow velocity; the enhanced tube did improve the heat 
transfer, however, the enhancement effect was weak in 
the case of high x, 


In laminar flow region, h decreased along the flow 
direction and reached a nearly constant value in the fully 
melted region; Nujocq increased with x, and was almost 
independent of heating power; In turbulent flow region, 
h decreased along the flow direction in the supercooling 
region, and then decreased in the melting region, then 
increased in the fully melted region; higher heating 
power caused higher h in the last region; Num of 5 wt% 
was higher than that of 10 wt% MPCS at same Re 


hm increased with vs; hm increased with x, when vs was 
0.05 kg min~', but decreased with x, when vs was high 


Based on the dimensionless wall temperature, MPCS had 
better heat transfer performance than that of water 
except at very small dimensionless axial distance 


Nutocai decreased along the flow direction; Nujocai was 
close to that of water; the heat transfer was enhanced by 
TiO} nanoparticles, and the enhancement ratio increased 
with x, 
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h decreased significantly at the entrance of the heating 


tube, thereafter decreased smoothly; h was 


1.82m 20-50kgh~! 


Copper tube, D;=10 mm, L 


dp=1-20 um, 


Paraffin, polymer 


Delgado 


10 wt% 


Xs= 


et al. [156] 


approximately 25% higher than that of water at the same 


Vs 


Compared to water, MPCS was beneficial since which 


could achieve lower wall temperature 


0.1-0.4 kg min7! 


Nine parallel minichannels, 


2 mm x 4.2 mm x 350 mm 


dp=1.5-12 um, 


n-Eicosane; 


Dammel 
and 


polymethylmethacrylate x,=10-20 wt% 


Stephan 
[125] 


10m 0.8-1.1 ms™! 


h increased in the melting region, and then decreased in 
the fully melted region; h increased with Re; h increased 
slightly with x, at the same Re, but decreased distinctly 


with x, at the same vs 


Copper tube, D;=10.5 mm, L 


7, 


dp=2-10 um, Xs 


11 wt% 


Octadecane, gelatin 


Tumuluri 


et al. [15] 


h of the blend was always higher than that of MPCS 


MPCS with multi-walled 11 wt% MPCM; 


carbon nanotubes 


1.1 wt% MWCNT, 
d,p=10-100 nm, 
Lp=0.5-40 nm 
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in Fig. 27(c). There were ice fraction thresholds at about 10- 
15 wt%, after which the heat transfer coefficient increased dras- 
tically. The authors believed that the transition of ice slurry flow 
pattern from homogeneous to heterogeneous caused the accu- 
mulating of ice particles near the top tube wall and the direct 
contact of ice particles with tube wall enhanced the heat transfer. 
On the other side, Illán and Viedma [133] observed that the 
Nusselt number decreased with the ice fraction at the same 
Reynolds number, as shown in Fig. 27(d), and same conclusions 
were reached when using corrugated tube. 

The results given by Yamagishi et al. [9] showed hoca of MPCS 
flowing through heated tube was higher than that of pure water 
and increased with the slurry fraction, however, the comparison 
should be conducted at the same Reynolds number. The flow 
turbulence of slurry would be damped as increasing particle 
fraction, which caused the reduction of heat transfer coefficient, 
therefore hiocaı of MPCS decreased as increasing particle fraction 
under the condition of same flow velocity. Ma et al. [28] noticed 
the same phenomena when studying TBAB CHS flowing through 
the heated tubes. However, other researchers reported different 
results about the heat transfer of MPCS in the heated tubes: Wang 
et al. [12,13] found the Nusselt number of 5 wt% MPCS was higher 
than that of 10 wt% MPCS even at the same Reynolds number 
(2100-3452); Rao et al. [152] reported the mean heat transfer 
coefficients of MPCS in a minichannel, which increased as slurry 
fraction at low flow velocity (0.05 kg min~') but decreased in the 
case of high flow velocity (0.15-0.35 kg min~'); Tumuluri et al. 
[15] measured the heat transfer coefficient of 7 wt% and 11 wt% 
MPCS in heated tube, showing hoca of 11 wt% MPCS was nearly 
the same as that of 7 wt% MPCS if the Reynolds number was the 
same. However, hiocaı decreased distinctly with slurry fraction at 
the same flow velocity. 

Song et al. [27] obtained quite similar variation of heat transfer 
coefficient of type A TBAB CHS flowing in heated tube to that of 
ice slurry reported by Niezgoda-Zelasko [145], as comparing 
Fig. 28(a) to 27(b). At a certain flow velocity and as the increase 
of slurry fraction, the flow of TBAB CHS was gradually damped 
from turbulent state to laminar state, therefore the heat transfer 
coefficient changed little with low slurry fraction (x, < 10 wt%) 
since the flow velocity dominated the heat transfer performance 
in the turbulent flow region. The re-laminarization occurred at 
10 wt% < xs < 15 wt% which caused large reduction of heat trans- 
fer coefficient, thereafter the heat transfer coefficient increased in 
the laminar flow region. Kumano et al. [31] also reported the 
similar results of TBAB CHS, as shown in Fig. 28(b1) and (b2), 
however, the re-laminarization occurred at the slurry fraction of 
about 15-17 wt%. 


3.4.3. Influence of flow passage geometry on heat transfer 

Niezgoda-Zelasko and Zelasko [148] discussed the impacts of 
tube diameter on the heat transfer coefficients based on the test 
results of ice slurry flowing through circular tubes and rectan- 
gular channels. Smaller hydraulic diameter enhanced the micro- 
convection effect of ice particles and therefore improved the heat 
transfer coefficient. However, other factors bring disadvantages 
for the heat transfer in the case of smaller tube: Reynolds number 
at the same flow velocity was lower than that in the case of larger 
tube, and the melted near-wall layer was thicker which hindered 
the heat transfer between tube wall and the fluid core. Thus, the 
authors believed that the smaller hydraulic diameter the smaller 
increase of heat transfer coefficients compared to the single- 
phase carrier fluid. Results reported by Kumano et al. [150] about 
ice slurry also supported this conclusion. 

The application of corrugated tube by Bédécarrats et al. 
[123] enhanced the heat transfer of the flowing ice slurry, about 
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Table 7 
Heat transfer of CHS. 
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Authors Initial Slurries Geometries Flow rates Correlations Comments 
solution 
concentration 
Takao TBAB CHS Plate-fin heat 1-11 kg min7! hm generally increased with vs; hm was 
et al. with cold exchanger, larger than that of water at 1-6 kg min7! 
[21] carry capacity Copper tube, flow rate, however, was smaller at higher 
of 64- D;=10 mm flow rate; hm was little smaller when the 
95 kJ kg! cold carry capacity was higher 
Song 30 wt% Type A TBAB Stainless steel 0.35-6.8ms—! Laminar flow: hm generally increased with vs; hm 
et al. CHS, x,=0- tube, D;=14 mm, Num = 1.231 x 107~>Re?-8606 py0.7073 | increases slightly or was near constant in 
[27] 18.4 wt% L=2.8m Turbulent flow: the turbulent flow region; as xs increasing 
Num = 5.254 x 10~4Re%9097 pr1.1202 to a certain value, hm dropped largely, and 
the flow was re-laminarized; hm increased 
with x, in laminar flow region 
17.3 wt% Type B TBAB Laminar flow: Variation of hm was similar to type A CHS, 
CHS, xs=0- Num = 1.798 x 107? Re®-6648 py1-3309 except that hm in turbulent region 
26.2 wt% Turbulent flow: increased rapidly with x, 
Num = 6.085 x 10~4Re®9130 pr1.2061 
Ma etal. 10 wt% Type B TBAB Copper tube, 0-2.5ms~'in Laminar flow: h was enhanced compared to water at the 
[28] CHS, ps=0- D;=6, 14mm, 6 mm pipe, 0- Nu(x) = 0.2368(Reyp Pr)*°873 Ste-0.0795 same Reynolds number; h decreased at the 
20 vol% L=1.1m 1.2ms™! in Di 02645 fd y 00182 melting region and increased at the 
14 mm tube (2) @) following solution region; h increased with 
Transitional and turbulent flow: Ps and Reynolds number, while ps was 
0.7093 ¢,. 0.0498 more important for laminar flow and 
Nu(x) = 0.0219(Remp Pr) ete Reynolds number dominated the heat 
(2) (2) transfer in turbulent flow; the influence of 
x D; heating power was distinct in the case of 
6 mm tube 
Maand 10wt% Type B TBAB Plate heat 2.5-13 L min“! Nu, = 0.338Re9.657 pr1/3 (1 4. Awa 0.037 Migcai fluctuated before the transition from 
Zhang CHS: y;=0- exchanger, 10 local = MR ( bs seat) slurry to aqueous solution, and increased 
[33] 17.5 vol% plates, D,=4 mm, sharply when turning into solution; Noverau 
W=72 mm, basically increased with the inlet water 
H=310 mm temperature, however, it dropped as the 
increase of p; until an un-fully melting 
happened 
Kumano 22.5 wt% Type A TBAB Stainless steel 1000 < Re < 7500 Type A TBAB CHS: In laminar flow region, h increased with xs, 
et al. CHS, 0- pipe, Dj=7.5, Nu = 3.33Gz9;18x942 (D/l) oaa but decreased along the heating tube, 
[31] 20 wt% 10.2 mm, L=1 m Type B TBAB CHS: especially at the tube inlet, and the effect 
8.5 wt% Type B 0.32 of Re was not significant, the ratio of 
Hat CHS, Nu = 6.35GZM x57? (D/) measured Nu to fié calculated single- 
0-20 wt% phase Nu was about 1-2.5 and increased 
with xs; in turbulent flow region, h was 
almost constant with x,, but dropped at 
about 15 wt%, then increased slightly at 
high xs, h changed little along the heating 
tube and increased with Re distinctly, the 
ratio of Nu was almost 1 with all x, 
Suziki 35 wt% 12 wt% Copper tube, 0.5-1.5 m s7! With the same Re, Nu of slurry was higher 
et al. ammonium ammonium i=13 mm, than that of solution and water; Nu 
[157] alum alum hydrate L=2m decreased as the addition of surfactant 


slurry 


caused by the laminarization 


1.5-2.5 times heat transfer coefficient was achieved in compar- 
ison with smooth tube. The same conclusion was reached by 
Alvarado et al. [11] about MPCS. They compared the heat transfer 
coefficient of MPCS flowing through both regular and enhanced 
8 mm tubes at the same flow velocity, and found that the average 
heat transfer enhancement was about 30% with 7 wt% slurry. 
However, the enhancement reduced when the slurry fraction 
increased due to the turbulence damping. 


3.4.4. Heat transfer characteristics of PCS in heat exchangers 

In addition to in the heated tubes, heat transfer performances 
of PCS in the heat exchangers have also been investigated. 
Fig. 29(a)-(d) shows the typical heat transfer coefficients of ice 
slurry in the plate heat exchangers and double-tube heat exchan- 
gers, all with hot water in the other side of heat exchangers. As 


the increase of ice fraction from 5 to 22 wt%, Bellas et al. [135] 
obtained nearly constant values of the overall heat transfer 
coefficient in a plate heat exchanger, as shown in Fig. 29(a). 
Moreover, the overall heat transfer coefficient reasonably 
increased with the flow rate. The authors believed that the overall 
heat transfer coefficient was the function of ice crystal size, 
Reynolds number and the flow passage geometry. Different from 
the results reported by Bellas et al. [135], the overall heat transfer 
coefficient in plate heat exchanger measured by Norgaard et al. 
[136] increased with the inlet ice fraction, as shown in Fig. 29(b). 
As shown in Fig. 29(c), Lee et al. [146] found that the heat transfer 
coefficient of ice slurry flowing through double-tube heat exchan- 
ger showed weak dependence on ice fraction when the flow rate 
was low, and there was a sharp increase when the ice fraction was 
higher than about 10%. However, when the flow rate was high, the 
heat transfer coefficient increased gradually as the ice fraction 
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Fig. 23. Ratio of Nusselt number of ice slurry to that of the chilled water [112]. The references in the figure are those in [112]. 
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increasing from 0% to 26%. The authors also pointed that the 
impact of hot water temperature on the heat transfer coefficient 
could be neglected in spite of the increase of heat transfer rate 
with hot water temperature. Bédécarrats et al. [123] presented 
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Fig. 25. Local heat transfer coefficients of MPCS flowing turbulently through 
heated tubes under constant heat flux, (a) Yamagishi et al. [9], 10.1 mm tube and 
(b) Wang et al. [12], 4mm tube. 


the variation of heat transfer coefficients of ice slurry flowing 
through double-tube heat exchanger with smooth and corrugated 
inner tubes, and the typical results in smooth tube is shown in 
Fig. 29(d). It was noticed that the variation of heat transfer 
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coefficient was similar to that of the pressure drop shown in 
Fig. 20(b). At lower flow rate range, the flow of ice slurry with 0% 
ice fraction was in the laminar region, therefore the increase of ice 
fraction led to more ice particles melting near the heat transfer 
wall and hence the heat transfer coefficient increased; at higher 
flow rate range, the flow of ice slurry with 0% ice fraction was in 
the turbulent region, and the turbulence decreased with ice 
fraction until the flow changed to laminar flow, therefore it was 
noticed that the heat transfer coefficient increased with ice 
fraction in the turbulent region and followed by a reduction in 
the transition region, and then a more rapid increase in the 
laminar region, as shown in the figure. The authors attributed 
the rapid increase of heat transfer coefficient in laminar region to 
the thinner ice-free thermal boundary layer which was caused by 
the enhanced particle-wall interactions with higher ice fraction. 
Higher hot water temperature was found to result in higher heat 
transfer coefficient both in laminar and turbulent regions except 
in the transitional region. Fernández-Seara et al. [118] investi- 
gated the heat transfer performance of ice slurry in a terminal 
fan-coil unit. Results showed the heat transfer was enhanced as 
using ice slurry instead of chilled water by the factor of 3.7-4.6, 
which depended on the fan rotation rate as well as the ice 
fraction. The dominated thermal resistance of the overall heat 
transfer was on the air side, where the ratio was in the range of 
80-88%. However, Illán and Viedma [158] did not recommend the 
direct application of ice slurry in fan-coil unit, since the frost 
formation, the partial blockage of fan coil inlet collector and the 
excessive diameter ratio (d,/D;) displayed the disadvantages of ice 
slurry compared to the single-phase flow. 

Ma and Zhang [33] investigated the heat transfer of TBAB 
CHS in a plate heat exchanger, and concluded that the local heat 


transfer coefficient was highly dependent on the flow velocity 
but slightly dependent on the slurry fraction, and hiocaı would 
increased drastically after the slurry turning to solution due to 
the low viscosity and high Re of solution. Meanwhile, the 
overall heat transfer coefficient decreased with the slurry 
fraction both at a certain flow velocity and inlet hot water 
temperature. 


3.4.5. Heat transfer correlations 

Except Knodel et al. [132] who developed the heat transfer 
correlation of ice slurry flowing through heated tube simply by 
multiplying a factor on Petukhov correlation of single-phase fluid, 
other researchers considered the heat transfer coefficient or 
Nusselt number of ice slurry flowing through heated tube or heat 
exchanger as the complicated function of Reynolds number, 
Prandtl number, ice fraction, ratio of the particle diameter to the 
size of flow passage and the ratio of the viscosity of bulk fluid to 
that of the fluid at the wall and so on. For the correlation of local 
Nusselt number, Graetz number (Gz=RePrD)/x) is widely utilized 
[112,147-149]. The index of Gz in the correlation is positive, which 
indicates the increment of Nusselt number with Reynolds number, 
Prandtl number as well as the entrance distance. The index of ice 
fraction is also positive to reflect the impact of ice particles on heat 
transfer. As aforementioned, the influence of tube diameter on heat 
transfer lies in both the micro-convection effect of ice particles and 
the near-wall thermal layer, therefore the index of diameter ratio 
(d,/D;) might be positive and negative [133,145,148,150] under 
different situations. Niezgoda-Zelasko [145] and Niezgoda-Zelasko 
and Zelasko [148] recommended a term of ((AH /CpAT)(Ax; /100)) 
to reflect the influence of ice melting on the heat transfer when 
developing the heat transfer correlation, which showed good 
agreement with experimental results. 

Similar to ((AH/C,AT)(Ax;/100)), Charunyakorn et al. [159] 
proposed the utilization of Stefan number, (Ste= (Cpq,,D;/ 
(2/s))/(XsAH)), to reflect the influences of latent heat and heating 
power on the heat transfer in developing the heat transfer 
correlation of MPCS flowing through heated tube, and Goel et al. 
[151] validated the importance of Ste on the heat transfer of MPCS 
experimentally, therefore Wang et al. [13] used this parameter to 
develop the heat transfer correlation of MPCS in heated tube and 
the relative error was shown to be within + 10%. Ma et al. [28] 
also utilized this term for TBAB CHS flowing through heated tube 
and the similar term of (1 + Ax;AH/(C,AT)) for TBAB CHS flowing 
through plate heat exchanger [33]. Good results were achieved 
using these latent heat terms. 


4. Application of PCS to refrigeration and air conditioning 
and performance evaluation 


4.1. Ice slurry 


Because of the numerous advantages, ice slurry is undoubtedly 
promising and can be used in various applications, including air 
conditioning, refrigeration, process cooling and so on. Many ice 
slurry systems were installed in Japan, which mainly concen- 
trated in air conditioning system with thermal energy storage, 
while the applications in Europe and USA cover all the fields. The 
typical utilization of ice slurry can be divided into the following 
ways: direct contact cooling application, typically represented by 
fish cooling system; direct application, ice slurry flows directly to 
the load side and exchanges the thermal energy with the high 
temperature air or water; direct carrier liquid application, only 
liquid phase in the slurry is pumped to the load side to avoid the 
blockage in the narrow flow passages; indirect application, ice 
slurry is only employed as the thermal storage material while the 
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tube, ø is the inlet ice fraction [133]. 


cold energy is first exchanged and delivered to the load side by 
another coolant, such as chilled water. 

A typical application of direct ice slurry application system in 
large-scale building is the CAPCOM building in Osaka, Japan 
[3,160]. The building is 20 stores and totally has 16,784 m° floor 
area, where the 2nd to 16th floor are cooled by ice slurry. The ice 
slurry is generated and stored at the 17th floor, as shown in 
Fig. 30. There are two ice slurry units driven by gas driven 
absorption refrigeration system with 272 kW cooling capacity 
each. The ice slurry is generated by a shell-and-tube ice generator, 
and ice is formed inside the tube and is then stored in 96 m? 
storage tank. Each unit has three alternating operated heat 
exchangers, and two for freezing and one for defrosting. 8 mixers, 
which are totally 55 kW, are installed in the tank to agitate the 
mixture of ice and water, then 40 wt% ice slurry at the tempera- 
ture range of —2.4 to —2.6°C is obtained. The ice slurry is 
pumped from the tank first through an ice slurry separator to 
adjust the slurry fraction to 20 wt%, and then reaches the air 
handing units on the floors of 2-16. Due to the low temperature 
and large cooling capacity of ice slurry, the air distribution 
temperature in the handing unit is designed to 12 °C which is 
lower than the normally-used temperature at 15 °C, therefore the 
air flow rate is reduced by about 22% which leads to capital and 
operation cost reduction. Compared to the predicted performance 


of an equivalent conventional ice thermal storage system which 
has higher air distribution temperature, the total energy con- 
sumption of ice slurry system during the summer of 1996 
achieved 4% reduction resulted from the energy saving in 
hydraulic and air distribution. 

An indirect ice slurry system was built for a 190,000-square- 
foot basketball arena and athletic complex cooling usage at 
Virginia Commonwealth University in USA [161], which has a 
seating capacity of 7500 people. The system schematic is shown 
in Fig. 31. The base aqueous solution used for ice slurry genera- 
tion is 7% glycol solution, of which the freezing temperature is 
—2.2 °C. As shown in the figure, the solution is pumped to the 
ice slurry generator which has the capacity of 380t (about 
1336 kW), and —8.3°C primary refrigerant is utilized for ice 
generation. The ice slurry is generated and stored at —3.3 °C. 
However, the ice slurry is not pumped to the load side directly, 
only the chilled solution at 1.1 °C is extracted from the tank to 
the building heat exchanger to cool water, which is located 
about 100 yards from the ice slurry system. After heat exchan- 
ging, warm solution at 11.6 °C returns. Thereafter, 2.2 °C water is 
supplied to the air handing unit to satisfy the proper fresh air in 
the 7500 seat arena. Compared to the conventional 5-7 °C 
coolant, the supply of 2.2 °C water causes the scale reduction 
of many devices, such as air coil, fan motor, water pump, duct 
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tube, Kumano et al. [31]. 


and pipe etc. The annual operation cost is reduced by approxi- 
mately $75,000. 

Sunwell Technologies Inc. [162] has supplied a number of direct 
contact ice slurry systems for fish processing and preservation as 
well as for supermarket display usage. The named Deepchill 
Variable-State Ice System can provide 0-100 wt% ice slurry with 
0.25-0.5 mm round crystals at the temperature range of —30 to 
0°C, which is easily pumpable. All components that employed in 
the ice slurry generation, storage and delivery are constructed by 
medical grade stainless steel or food grade plastic composites, 
while the slurry is always sealed and never exposed to outside 
contaminants. Hence, the slurry is considered as clean and reliable 
enough for direct contact with food. 25-40 wt% ice slurry is usually 
utilized and can cool the target goods to —2 to 0°C in several 
minutes. 

Ice slurry system constructed in a large European institutional 
kitchen in France [4,163] provides a complete cooling chain in 
this kitchen, including refrigerated room, refrigeration system, 
refrigerators, recooling vessels, the food service conveyor belts, 
plate dispensers and cooling tray transport trolleys. This system 
employs 16 scraped surface ice slurry generator of 185 kW total 
cooling capacity and 3 cylindrical storage tank of 22 m? to 
generate and store ethanol/water ice slurry. It was found that 
the usage of ice slurry could provide a constant low temperature 
for food processing and satisfy all the cooling requirements in the 
large kitchen 


4.2. MPCS 


An MPCS thermal energy storage system was constructed for 
cooling in Narita Airport in Japan [18] without replacing the 
existed 2500 RT (about 8792.5 kW) electrical turbo chiller. This 
project aimed to improve the decrease in capacity and perfor- 
mance caused by refrigerant substituting for environmental con- 
servation. Non-toxic and phase change temperature controllable 


material, n-Paraffin waxes, was selected as the PCM. The slurry is 
the mixture of water and 2 um microcapsules, and the melting 
temperature and the corresponding latent heat of which are 8 °C 
and 75.9kJkg~', respectively. The storage tank is 24.7 m in 
height and 7.4 m in diameter. Cold energy is stored during night 
time. MPCS at 11.5 °C is pumped by ordinary centrifugal pumps 
from the top of the storage tank, and is cooled in a plate-fin type 
heat exchanger by 3.5 °C coolant, as shown in Fig. 32(a). Phase 
change process occurs, the slurry is then cooled to 4 °C and stored 
at the bottom of the tank. The stored cold energy is released 
during day time as shown in Fig. 32(b). By the heat exchange with 
cold MPCS, the coolant of 5.1 °C is obtained and pumped to the 
load side. COP of the chiller and the system during the cold 
storage process is about 5.4 and 3.2, respectively, which is lower 
than the COP of 5.9 and 4.4 when there is no energy storage. 
Compared to the old storage system using external melting ice, 
though the thermal storage density of the used MPCS system is 
about 60% smaller, the operation cost is 32% lower. 

Griffiths and Eames [16] investigated the performance of MPCS 
as the heat transfer medium in a test chamber containing a 
chilled ceiling panel, and the schematic diagram of which is 
shown in Fig. 33. The authors employed the MPCS with a melting 
temperature of 18°C and 2-8 um microcapsules produced by 
BASF, and the phase change temperature is close to the working 
temperature of chilled ceiling at 16-18 °C. The authors compared 
the performance of chilled water and MPCS with 40% concentra- 
tion. To maintain the test chamber at a constant temperature of 
19 °C, flow rate of 0.7 Ls~! was required when water was utilized 
as the heat transfer medium, while only 0.25 Ls~! was required 
for 40% MPCS to maintain the space temperature at 20-21 °C, 
since the average specific heat of the slurry in the temperature of 
16-17.5 °C was as high as 9 kJ kg~' K~', with sometimes a value 
of 11-12 kJ kg~' K~'. Through the continuous operation for over 
four months, the MPCS with 40% concentration was validated as 
heat transfer medium on chilled ceiling panel instead of water, 
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and there was neither deposition of microcapsules in pipe nor the 
degradation of the slurry, and the system pumping power was 
reduced due to small flow rate. However, the authors noticed 
that the control procedure for cooling the slurry should be 
optimized. 


4.3. TBAB CHS 


Takao et al. [21] introduced the operation results of the TBAB 
CHS air conditioning system in the office building of NKK Keihin 
Works, which was constructed by the JFE Engineering Corpora- 
tion. Total area of the two floors building is about 4000 m2, where 
1700 m? is the cooling object of TBAB CHS. As shown in the 
system sketch in Fig. 34, an absorption refrigerator is utilized. The 
cold carry capacity of the employed TBAB CHS is approximate 
2 times of that of chilled water, which is about 59 MJ m~’, 
therefore the coolant flow rate was reduced by 50%. The heat 
exchange between TBAB CHS and air inside the room was under- 
taken by plate fin heat exchanger, and the diameter of the cooper 
tube for TBAB CHS was reduced to 10 mm because of the low flow 
rate. The overall heat transfer coefficient in the heat exchanger 
was about 250-300 W m~? K~! as the CHS flow rate varying from 
1 to 10 kg min~!, which was higher than that of water when the 


flow rate was smaller than about 6.0 kg min~'. However, the 
overall heat transfer coefficient was lower at higher flow rate 
range. Compared to the utilization of chilled water, though there 
was a 20% increase of the pressure drop of TBAB CHS, the 
operation pumping power measured in the representative day 
of August of 2000 was reduced by 47% due to the lower flow rate, 
while the total energy consumption had a reduction by 21%. 

Ogi [164] reported the thermal energy storage air conditioning 
system using TBAB CHS in Kawasaki Underground Azalea which 
has a total area of about 56,000 m?. The project was also 
conducted by JFE Engineering Corporation and was finished in 
2008. The system schematic diagram is shown in Fig. 35. The 
system was modified based on an existing air conditioning 
system, which consisted an air handling unit (AHU) system using 
3 heat pump of 550 Rt (about 1934 kW) each, and a fan coil unit 
(FCU) system with a centrifugal chiller of 500 Rt (about 1758 kW). 
Thermal energy storage by water was employed by the FCU 
system using a storage tank of 700 m”. When TBAB CHS was 
introduced to the FCU system instead of water, the centrifugal 
chiller was substituted by a screw chiller with 500 Rt cooling 
capacity and the COP was improved from 4.19 to 4.44. Since the 
cold carry capacity of TBAB CHS was about 2 times as large as that 
of the chilled water, the stored cold energy was improved from 
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Fig. 31. Schematic diagram of indirect ice slurry application in university arena 
[161]. 


840 Rth (about 2954kWh) by water to 1650Rth (about 
5803 kWh) by TBAB CHS without changing the storage tank. 
Compared to the original water storage system, the new TBAB 
CHS storage system led to a 40.17% reduction of electricity 
consumption and a 40.28% reduction of CO, emission. 

Moreover, all the other TBAB CHS projects conducted by JFE 
Engineering Corporation, including the systems of Takenaka 
Corporation (2004), Tsurumi office building of JFE Engineering 
Corporation (2005), Kuraray Kurashiki Office Production Technol- 
ogy Center (2007), NTN Okayama Plants (2009), achieved more 
than 20% energy saving [165,166]. 
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Fig. 32. Operation strategies of the MPCS system during (a) night time and (b) day 
time [18]. 


5. Conclusions 


In this paper, the application of phase change material slurries, 
i.e., ice slurry, MPCS as well as CHS in the secondary loop 
refrigeration and air conditioning system have been reviewed, 
and the utilizations of PCS are expected to facilitate the energy 
saving and environmental conservation. The paper mainly focuses 
on the generation method and storage strategies, the flow and 
heat transfer characteristics and the applications of these slurries 
in real projects. There are many common generation methods of 
ice slurry in which the necessary technologies are used to prevent 
the generated solid crystals from adhering to the heat transfer 
surface to deteriorate the heat transfer, e.g., using scraper or rod 
or fluidized solid particles to detach the adherence, using the 
supercooling method to avoid the crystals forming on heat 
transfer surface, using direct contact heat exchange to improve 
the heat transfer, and so on. The producing process for MPCS is 
relatively complicated and higher cost, which may limit its 
application. The homogeneous storage of ice slurry and CHS 
benefits the extraction but has limitation on storage density, 
while the heterogeneous storage can achieve high storage density 
but requires additional power for extraction and the crystals size 
may become larger. Though PCS has been widely investigated and 
utilized, the fundamental investigations on flow and heat transfer 
characteristics are still divergent because of the complex influ- 
ences by many factors, therefore further studies are still required. 
Moreover, the application of these slurries show quite positive 
significance on energy saving and environment conservation. 
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Fig. 33. Schematic diagram of chilled ceiling in the testing chamber [16]. 
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Fig. 34. Schematic diagram of air conditioning system using TBAB CHS in the office building of NKK Keihin Works [21]. 
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Fig. 35. Schematic diagram of air conditioning system using TBAB CHS in Kawasaki Underground Azalea [164]. 
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However, MPCS and CHS are still not thoroughly investigated and 
applied, which suggests that more fundamental researches are 
still necessary for MPCS and CHS. 
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